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STABILITY AND STABILITY SWITCHING IN 
3-SPECIES MODEL WITH DELAYS 

IDRISH H. ELMABRUK* AND H. L. ARORA* 

(Received 02.01.2002) 

ABSTRACT 

Sufficient conditions for the stability of the positive equilibrium of a 3-species Lotka- 
Volterra-type system with density-dependent death rates modeling the interactions of two 
competing predators living exclusively on a common prey with delays, are obtained by 
constructing a suitable lyapunov function. An estimation of delay length for which equilibrium 
stability holds is obtained. We show the existence of some critical value £ of the delay a in the 
predation and biomass conversion terms at which the stability switches to instability and the 
system remains unstable for & > q . 

Key words: Predator-Prey System, Equilibrium, Lyapunov Function, Delay, 
Stability-Switching. 


INTRODUCTION 

One of the most interesting questions in ecological models concerns with local 
stability of the positive equilibrium of the model. The question becomes even more 
interesting when we consider time delays in the growth dynamics of the interacting 
populations of various species. Delays can arise from a great variety of causes and 
undoubtedly always present to some extent. In order to, realistically, model a real life 
system of inter-specific (and intra-specific) interations, one generally needs to take 
into account the influence due to the past states of such a system, or the waiting 
period resulting because of physiological processes. These usually mean that a real 
system should be modeled by differential with time delays. Indeed, such a view is 
gaining popularity across the board of science and engineering, as evidenced by the 
rapid increase in the number of mathematical models of delay differential equations 
studied in various scientific disciplines, e.g., see the numerous references cited in [1 ]. 
This is specially true in the area of mathematical biology, see [4,8,11]. 


* Department of Mathematics, University of Garyounis, Benghazi, Libya. 
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Considering the historical Lotka-Volterra 2-dim system without delays, that is 



— = -dH + cHP 
dt 


( 1 . 1 ) 


where a, b, c, d are positive constants, it is very well known that (1.1) is unrealistic for 
various reasons. We, for simplicity, nondimensionalize( 1.1) by 

, x cH(t ) . . bP(t) d 

putting m(t) =—-—,v(t) =-, x = at,a = —, then (1.1) becomes: 

da a 


du 

dx 


= «( 1 - v) 


rfv 

dx 


= av(i/-l) 


( 1 . 2 ) 


The trajectories of (1.2) are given by 

au + v-vln(u a } = N,whereA^ > 1 + a, (1.3) 

and in the phase plane are closed (periodic solutions) as shown in Fig. (1.1), a typical 
periodic solution (w(t), v(t)) is shown in Fig. (1.2). A major inadequacy 



Fig. (1-2) 

of (1.1) is that its solutions are not structurally stable. Suppose for example m( 0) and 
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v(0) are such that u and v for t > 0 are on the trajectory N 4 which passes close to the 
u and v axes. Then any small perturbation will move the solution onto another trajectory 
which does not lie everywhere close to the original one N A . Thus a small perturbation 
can have a very marked effect, at the very least on the amplitude of the oscillation. In 
other words, this high sensitivity to the smallest change in the model leads to a model 
with essentially different species behaviour. For instance the positive equilibrium point 
of the system which is a center changes, under small pertubations, to a focus. For a 3- 
species Lotka-Volterra-type model without delays, as show in [5], there exist Hopf- 
bifurcating periodic orbits, but no conclusion about their orbital stability can be drawn 
by using the available procedures given in [10,13] and others. 

Thus, instead of relying on the instantanous density of one species (prey, say), 
because young preys do not participate in the reproductive process prior to sexual 
maturity, together with the fact that the predators are not usually successful in hunting 
young preys which do not wander away from their safe havens, we can get a more 
realistic model if we take into account density of one species (prey say) in the past, 
that is a system with delays. 

There are several examples of 3-species communities such as linear chain of 3- 
species, see Fig. (1.3) (i); two preys and one predator, see Fig. (1.3) (ii) and two 
predators interacting with one prey as in Fig. (1.3) (iii). 



Fig. (1.3) 

Since the question of two predatory populations consuming a single prey 
population has occupied an important place in ecological literature, thus in this paper 
we have a differential system modeling the interaction of two predators feeding upon 
a common prey. We have incorporated discrete delays t > 0 in the logistic maturation 
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term of the prey species and a > 0 in the predation terms and biomass conversion 
terms. Both predatury populations have been considered with density-dependent death 
rates. By using Lyapunov function technique as discussed for 2-dim predator-prey 
system in [9], we have obtained an estimation of delay's length for which stability is 
preserved, we have shown the existence of a critical value a of o (with t = 0) at 
which stability of the unique positive equilibrium switches to instability and the system 
remains unstable for o > u • 

MODEL WITHOUT DELAYS AND LOCAL ASYMPTOTIC STABILITY 

An intermediate type 3-species system, see [7], modeling the inter-species 
interaction of two predators and one prey is of the form: 

y = y[-g(y)+c<h{x)] 

z = z[-h(z) + dq 2 (x)] 

with 

x(0) = x 0 >0 , y (0) = y 0 > 0 , z(0) = z 0 >0 , 

J{x) is the specific growth rate of the prey x and generally is assumed to have the 

following properties: /[0,co)—> R\f eC'[0,oo),/'(x) < 0 ; there exists A'>0 such that 

J[k)= 0 where k is the environmental carrying capacity with respect to the prey x. 
q^x) and q 2 (x) are the functional responses of predators y and z repectively with 
respect to the prey x. 

Generally, it is assumed that q , ;[0,oo) -» —> [0,oo );q t eC'[0,oo);^(0) = 0,<?/(.v) > 0 

/= 1 , 2 . 

g(y) and h(z) are the density-dependent specific death rates of predators y and z 
respectively. It is assumed that g(0)>0, h(0)>0g'(y) > 0 for y > 0 and h '(z) > 0 for z >0. 

In this paper, however, we consider the folliwng special case of (2.1): 


( 2 . 1 ) 
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.V - .V 


V = V 


k ) 


a 


( 2 . 2 ) 


where the parameter A>0 is the above mentioned carrying capacity of the environment 
with respect to the prey .v which, in turn, has a logistic growth rate, (i.e. it does not 
grow unboundedly in the absence of predators), <=> 0 is the intrinsic growth rate of 
the prey .v in the absence of predation, a> 0 , b >0 are predation rates (i.e. functional 

responses of predators y and z respectively with respect to .v), —are biomass 
conversion rates. We, further; make the following assumptions: 

0 <u <b < 1 , k,e> b> a , 

and ka-eb>() if e* k ^2 3 ) 

and k+a>2b if e= k. 

The unique positive equilibrium E(x\y t ,z*) of (2.2) is given by: 


. k(€+a+b) , _ e(A -a) + k(b-a) 

€+2 k ' a(e+2k) 

, A(e-/?) + ka-b e 

/>(s+ 2k) 


and 


(2.4) 


It is straightforward to show that the eigenvalues of the variational matrix V(E) of the 
system ( 2 . 2 ) about £(x*,y*,z*) are given by the characteristic equation: 


f + P|X, 2 + + P x — 0 


(2.5) 
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where 


. e 

P\= x -+y +* 

f e\ e\ 

p, =xV 1 + - +xV 1 + -\ + yz* 
l k k 


* * * o 

p,=xyz 2 + 


Since P>0, P>0, and 


P,Pi - Pi -|(> + j )(*') V + *') + *'(l +J ][(/) ! + (z') ! 


**/* *\ 2 € • * « . 
+y z (y +z ) +—x y Z > 0, 

K 


then by Routh-Hurwitz criterion E[x ,y\z*) given by (2.4) is locally asymptotically 


stable. 


MODEL WITH DELAYS AND LOCAL ASYMPTOTIC STABILITY 

We incorporate discrete delays in the system (2.2) as follows: 


( \(t — \ 

x = x(t) € 1 --- - ay(t -ct)- bz{t -o) 

V k J 


y = y{t) -(1 + y( 0 ) + -x(t- ct) 

a 


z = z(t) -(1 -Z(t)) + -x(t-a) 
b 


with the initial conditions: 
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x(t) = 0 l (t),t e[-r, O],0, eC([-r,0],R*) 
y{t) = 0 2 {t\t e[-r,O],0 2 (t) e C([-r,0],/T) ( 
z{t) = 0 2 (t),t e[-rfi],0 2 (t) £C([-r,0],r) 


where r = max{T,cr}. 

Clearly the system (3.1) with (3.2) has a unique solution which exits and is positive for 
all t > 0.x denotes delay incorporated in the logistic growth term of the prey species. 
To make calculation simpler we have taken the delay in both predation tenns in the 
prey growth rate equation and in the biomass conversion terms in both predators 
growth rate equations to be the same (a >0 say), otherwise the analysis technique 
remains the same with different delays. 

Naturally E(x*'y’z % ) as given by (2.4) remains the unique positive equilibrium of the 

system (3.1). Let us introduce a new coordinate system through the following 
substitution: 

*,(/) = *(/)-/, y,(/) = y(r)-/ z t (t) = z(t)-z\ 

From(3.1), (x: 1 (t),y,(t),z,(r)) satisfies the system: 


M0=(M0+/)/ 

M0=(M0 +z >*( 


- 7 *1 (' - 0 - «.v, K (t - <0 

k 

I 1 

-y,(0 + -^(/-o) 

a 

-Z|(0 + jk(/-a)j 


(3.3) 


Remark 3.1. 

The stability analysis of E{xy , z *) is reduced to that of zero solution of (3.3). 
From (3.3), it follows that the variational system of (3.1) with respect to the positive 



64 


STABILITY AND STABILITY SWITCHING IN 3-SPECIES MODEL WITH DELAYS 


equilibrium E^x'y'z') is given by the linearized system of (3.3) as: 


*i (0 = **i 

i ; i (') = >'* 




1 


-y\(t)+-x i (t-o) 


% t \ 

Z\(t) = z' -*,(/)+-*,(/-o) 

o J 


(3.4) 


Now, we can rewrite (3.4) as follows: 


d_ 

dt 

d_ 

dt 

d_ 

dt 


v,(0 


— fx, (s)ds -a [ y^ds-b fz,(.yWs 
x k , J . , J _ , J _ 


h oo J 

• • 

y a 


i 

jxf(s)ds 




t-o i-a 


= y^ + -x l (t) 


= Z_ i (t) + -X,(t) 


x \{ t )~ a y\{t)h z \{t) (i) 

(H) 
(Hi) 


(35) 


Now, we construct a Lyapunov function w(t) so that the trivial solution (i.e. (0,0,0)) 
of (3.3) is stable and then we show solution (.xj(^),z 1 (/‘))->(0t 0,0) as / -> oo which 
implies that (0,0,0) is a locally asymptotically stable solution of (3.3), which 

consequently implies that E(x\y*,z') is locally asymptotically stable. 


Consider 

w(t) = X 

, w 1 (r) + /a 2v v 2 (/) + rVw 3 (/), 


(3.6) 

Where 

*»*. 

ii 

jE 

+ 

JS 

|s-» 


(3.7) 

For 

w i,W= 

1 l ^ j x, ( s)ds a J y x (s)ds - b j z, {s)ds 

X k i-x t-o i-o 

2 

(3.8) 
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T J ] x 2 \{l)dlds + a j \ y 2 (l)dlds + b \ jzf(l)dlds 

K I -T S /'OS /-os 

where w 2 (/) = w 21 (r) + w 22 (/) 



(3.9) 

(3.10) 


for 


for 

and where 

for 


W, 


=(') = 


MA+- 1 ^,( 5 )^ 

V a ,-o 


w a('H 7 + _ V) } \^{l)dlds 

> 2 d Jr-a s 


W »(0= W 3l(0 +W 32(0 


/ X I 2|(0 1 

^3.(0 = 


6 f -o 


(3.11) 

(3.12) 

(3.13) 

(3.14) 


*'«(') = 


(i + ±)jj, W 


(3.15) 


Remarks 3.2. 

Let X(t) = (.V|(t),^,(r),z,(/)). Thus by w(r) wen mean: w(t) = w(A')(r), and 

similarly for w 1 (/),w 2 (/),wj(/) f w ll (/),w 12 (/),w 2 l (/),w 22 (/),w jl (^) and w J2 (t). 

From the defination of w(t), it is clear that w(t) is positive definite. 

Next we compute the time derivative of w(t) with respect to (3.4). As shown in the 
Appendix we get: 


l M!l< 

+ «T (0 


* 


-2 + x 


2 e 

\T 


x + (a +6 )(t + ct)1 + — y*a(a + 2) + — c(b + 2)z* 


-2 + — f—(t + o) + 2a(tf +6)1 + — y* 
a \k J a 


+ b 2 z 2 (t) 


-2 + - 


: (t + a) + 2o(a +6) 


+ — z 
b 
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Thus 

(3.16) 

where it is assumed that: 



a = minj—a, a a 2 ,Zra 3 


(3.18) 


Thus (3.16) gives: 




(3.19) 


From (3.19) it follows that w(.x 1 (/),)> 1 (/),z 1 (r)) as defined by (3.6) is the sought after 
Lyapunov function. To prove local asymptotic stability of E[x\y\z*) it remains to 


show that: 




Now integrating (3.19) from T to t, we get 

+ + + w(T), t*.T 

T 

But (3.20) implies that x 2 (/)+y 2 (/) + z 2 (/) eZ^oo). 


(3.20) 
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We claim that }i™[**(0 + J’ 2 (0 + z ?(0] = 0 • 

Since w(t) is a decreasing function with strict minimum at (0,0,0), MO exists and 
it finite. 

This is because Jini[ -*| 2 (0 + T| 2 (0 + z i 2 (0] ~ ^ will according to (3.20) lead 

to a contradiction. Thus it follows that the (0,0,0) solution of (3.3) is locally 
asymtotically stable provided the following three conditions hold: 


2e 

{ k 


+ a+b k + 


x*(a+b) + y* — (a + 2) + z*—(6 + 2) 


a <2 


x 6 
V ak J 
*X ^ 


T + 


x ( e _ t 

— — + 2a + 2b H- 

a Kk ) a 


V bk J 


t + 


rV 


—h 2 a + 2b H— 

k J b 


a <2 

la <2 


(•) 

(«) 

(«) 


(3.21) 


Now, we determine values oft and a satisfying the three inequalities (3.2.1) (i), (ii) 
and (iii). 

|| * 

Since from (2.3) we have 7 < ~, and fromn (2.4) it follows v* > z and thus — >—, 

b a a b 

hence it is sufficient to find values t and a satisfying (3.2. l)(i) and (ii) only. Thus we 

have proved the following main theorem: 

Theorem (3.1) 

Assume that (2.3) holds. If the delays t and a satisfy the inequalities (3.2.1) 
(i) and (3.2.1) (ii), then the positive equilibrium E[x ,y* ,z*) is locally 
asymtotically stable. 

Let us consider the following illustrative example: 

Example 3.1. 

1 _ 5,6 . 19 . 4 

If we let e=k = \,a = -,b = ~, then from (2.4) we have * = — .T = —, 
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and z = - consequently the inequalities (3.2.1) (i), (ii) and (iii) become: 

O 


(2.67)t + (1.94)a<2 (/) 

(1.266)t + (5.176)ct<2 (//) ^ 

(l.055)t + (4.022)o < 2 (iii) 

Thus as long as 

0<t< .749 and] 

0 <a< .386, j 


(3.22) 


(3.23) 


the equilibrium E 


]9 4_ J_ 
24 *15*18. 


is locally asymtotically stable, see Fig. 3.1. 


Domain of attraction 


a 



Remark 3.3. 

It follows from (3.23)-that the delay a>0in the predation and biomass 
conversion terms is more destabilizing than the maturation delay t > 0 of the prey as 
far as this particular Lotka-Volterra predator-prey model is concerned. 

STABILITY SWITCHING 

In this section, we show the existence of some critical value a > 0 of the delay 
o > 0 (with t = 0) such that for a > a there occurs stability switching of the equilibrium 
£(x*,y*,z*) from stability to instability, in other words the equilibrium £(.v*,y*,z*) is 
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unstable for a > a. The discussion is based upon Theorem 1 of [2] and corollary 4.3. of [ 6 ]. 
We consider the system (3.1) with x = 0 and e — ^ = i, that is 

.v = -v(/)[(l - v(/)) - ay{t - a) - bz(t - a)] 

v= v(0 -(l+v(/)) + l.v(/-a) • (41) 

z = z{t) -(l + z(/)) + ^x(/-a) 

The characteristic equation for the linearization of (4.1) is 

A 3 + ak 2 + pA + y = (5 + AO)?- 2 '" (4.2) 

where, 


a = .v* + y‘ + z* 

p»,-(/ + ,-) + .vV 

* * * 

y - x y z 

5 = -2.y*vV 
0 = - .**(.)'* + Z*) 


(4.3) 


Let /■’(A.) — X + oX~ + pA + y 
(A(A) = -(A0 + 8 ) 

Since P and Q are polynomials with real coefficients, the condition (HI) of Theorem 1 
of [ 2 ] is satisfied. 

Now 

P(X) + Q(X) = A + oA , 2 + PA + y — A.0 — 5 = 0 or 

A , 3 +[x* +>’* + z*]A 2 +^ 2 .v*(/ +z*) + .vVj + 3.vV z * = 0 ^-4) 

Now (4.4) is same as (2.5) with (2.6). 

Thus as shown in section 2, £(.v\y\z’) is locally asymptotically stable. This implies 
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that P(X) + £?(>.) = 0 has no root in Re(X.) > 0. Thus condition (H2) of Theorem 1 
of [2] is satisfied. In order to utilize Theorem 1 of [2], let us define: 

= / +(ct 2 -20)/ +(p 2 - 2 oy -e 2 )/ +(y 2 -5 2 ) (4 ‘ 5) 

Where 


a 2 -2P = (/) +(/) +(/) >0 
F(0) = y 2 -5 2 = -3(//z’) <0 • 
0 2 - 2 ay - 0 2 =/z‘(/z*- 2 (/) 2 ) 


(4.6) 


Now Theorem 1 of [2] is stated as follows: 

Assume conditions (H 1) and (H2) hold. 

(a) If F(y) has no positive roots, then no stability switches occur. (This implies 
either the equilibrium point is stable for all <j > 0 or it is unstable for all a > 0 • 

(b) If F(y) has positive roots, stability switches may occur. At most a finite number 
of switches is possible. There exists 6 > 0 or such that the equilibrium point is 

unstable for all a > o . 

Now since a 2 - 20 > 0 and y 2 - 5 2 < 0, then by Descartes rule of sign applied to 

F(y) as given by (4.5), the possibility that F(y )>0 for ally > 0 (i.e. F(y) has no positive 
root) is completely ruled out. This implies that the stability of the equilibrium point 

£(/,/, z’) is not preserved for all cr>0- 
Thus we have following theorem: 


Theorem 4.1. 

Assume (2.3) (i.e. the equilibrium E[x ,/,z*) of (4.1) with <j = 0 is 

asymptotically stable). Then there exists a > o for which E(x ,y ,z ) is unstable 
for a > a . 

Remark 4.1. 

Prior ultimately shifting to instability for o > 6 , there may occur at most a finite 
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number of stability switches. 


Example 4.1. 

Consider the system (4.1). The associated linearized system has the characteristic 
equation: 


X 5 + aX 2 + pA + y - (0A + 8)e- Ua = 0. 


Let a = 



II * 

— ,e= a: = l, we have x 

O 


19 . 4 . 1 

24 ,y = \5 ,Z = 18' Thusfrom(4.3) 


we get 

a = 1.114, p = 02699, y = 0.0117, 5=-0.023457, 0 = -025499, (4.7) 

From (4.5), we have: 

F(y) = / + 0.70119/ - 0.018234/ - 0.0004107, 
which has positive real root y=0. 19724. 

This indicates that preservation of stability can not hold for all or > 0. 

Thus there exists a such that for all CT > 6 , the equilibrium is unstable. 


Remarks 4.2. 

The delay which could represent in our model influence of the past or waiting 
period in the proliferation of certain species due to (say) lag in some physiological 
process (because of natural or unnatural reasons), could have if longer than enough a 
destabilizing effect on the otherwise stable equilibrium coexistence of the ecosystem's 
species. 

Remark 4.3. 

In case the two delays x and a are kept in the model (4.1) as in case of (3.1), 
then the corresponding characteristic equation is: 

A, 3 +(/ +z’)A 2 +/z*A+^A 2 x* + Ajc*(A* +z*) + x*/z*Je‘ XT 

+[>/(/ +z*) + 2//z*Je” 2Xo = 0. 

However, discussion of the full analysis of stability switches for arbitary a and x 
remains an open problem, see [2,3], 

APPENDIX 


In this apprendix we prove (3.16). From (3.8), the time derivative of w,,(t) with 
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respect to (3.4) while using (3.5) is: 


dw u (t) 


-dt 


= 2 


x \( t ) ~ 6 \ 


x k 


J .y, (s)ds - a J.v, ( s)ds - b J r, (s)ds 




-2 €.vf(/) 2 a , . , . 2b , , , . 

—-^-r-v,(0 V, (/)- — .Y,(/)2,(/) 

A' .V .V A' 


+2 


j -V 1 (r).Y 1 (5)^ + ^|Y,(/);',(5)* + ^ j.v,(/)-,(*)* 


2 ea f 


J >-, (/)a, ( s)ds + 2a 1 J y x {t)y, (s)ds + 2 ab J y x (/f)z, (*)</.? 


2 eZ> f 


~ J-i(O v i(- y V v+2ab \-1 (0>',(• y V ' s+2t > 2 J-.(0-i( 5 K v 


-2 er, 2 (/) 2a , , , , 2b . . , > 
£ —-L_- r a, (/)y,(/)—r^,(0z,(0 

k A A A 


+ 


( e -v \ 


Y * A 


•vf(0v+ j.vf(i)^ .v, 2 (0<t+ 


€b 


,vr(/)a+ jz-’(.vVs + -^- >'f(/)t+ Ja 1 2 (5)c/. 


+a 


I I 

(0°+ JV (•*)<& +ab yf( 0°+ J *?(«)<& 


+- 


eb 


-f (/)t + |A, 2 (s)ds + a 6 z, 2 ( 0 a + J +b 2 z*{t)a + jz*{s)ds 


(A-1) 
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Define w,(/) = = w i 2 (^)(0’ (A-2) 

Where w l2 (/) is defined so that the integral terms in (A-l) get eliminated. Thus 




t I 

+/>J \z'(l)dlds 


(A-3) 


Using the fact that J x2 ( s ) d 

t-r 


d 

s + — 
dt 


t f 

| j x 2 (l)dlds 

t-r s 


- x '( t ) r , then from (A-l), (A-2) 


and (A-3), we get: 


fMO 

dt 


< -2ejj<) 2 a 


2b 


k x x x 




7 ) [*?(/)* + + 


+ y [*> (0 ct + z . 2 (*H+ y [y' M ct +*? W T ] 

+a 2 [y 2 {t)o + y 2 (r)a] + ab[y 2 (t)a + z 2 (/)a] 

+ y [ z ' M t + x ' W T ] + ab[z]{t)a + y 2 (t) ct] 
+b 2 [z 2 {t)a + z 2 (t)a] 


Thus, we have: 


dt 


zxlit)^ 

+y 2 {t)a 

+z 2 {t)b 


_4 + lfI + (l + i)( x + a ) 
x k 


—(x + a) + 2a(a +f>) 
k 

—(i + a) + 2 a(a + b) 
k 
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X X 


(A-4) 


Next, we define w v as in (3.11): 

Ti(') i 1 f 


V, ' 2 l(') = 


- L V^ + - J .v, (a)c/a- 
v a J _„ 


From (3.11) and(3.5)(ii) we have 


</w,|(r) 


dt 


= 21 


>.W, 1 f 


~+~ fx,(j)<fe 

>’ a ,-c J 


->i( 0 + -U (0 

a 


= +- -r *, (0>i (0 - - J y\ (0*. ( s ) ds +-A j ■*■ (O^i ( s ) ds 

y a y a L a j 


-2y;(t) 2 1 , w v 1 

<-^ + — TX t O.V. (0 + - 

y a y a 


i 

y;{t)o+ J*, 2 (j)ds 


<r 


i 

jr, 2 (/)a+ J jc, 2 (^)c/.v 


Define w,(/) = w 2 (X)(t) = w 21 (^)(t) + w 22 (A')(/) 

Where w„(/) is defined so that the integral terms in (A-5) get eliminated. Thus 

= j jx*(l)tllds 

V« « AV, 

From (A-5), (A-6) and (A-7), we get 


dw 2 {t) 


dt 


<Jcf( 0 CT 


1 2 

+ n 2 
a a 


+/(') 


2 a 
-r + ~ 

. y a 


+ ~\x t (t)y l {t). 
a y 


(A-5) 


(A-6) 


(A-7) 


(A-8) 


Similarly, we have: 

WjM = W 3(^X0 = W 3l(*X0 + W 3 2 ( *X0 


Where 
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and = (t + 7 t 1 I 1 x]{l)dlds, 

\0 O J,-ns 

and we get 




dt 


<xl(t)a 


1 2 
b + b 1 




-2 CT 

7" I T 


+ 7~W'M') 

ft z 


(A-9) 


From (3.6) and using (A-4), (A-8) and (A-9) we get (3.16). 
This completes the proof. 
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A TWO WAREHOUSE DETERMINISTIC ECONOMIC 
ORDER QUANTITY MODEL WITH FINITE RATE OF 
REPLENISHMENT AND STOCK DEPENDENT 
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ABSTRACT 

A deterministic inventory model is developed for a single item with two warehouses one 
owned by the management (abbreviated as OW) and other is a rented warehouse (abbreviated 
as RW) by considering stock dependent consumption rate. It is assumed that storage capacity of 
OW is limited and that of RW is Unlimited. The holding cost of RW is greater than OW and 
actual service of the customers is being done from OW, so it is economical to empty RW first. 
It is assumed that stocks are transferred from RW to OW according to K-Release Rule (KRR) in 
such a way that units in each transshipment is less than the capacity of OW. The model is set up 
by assuming finite replenishment rate and shortages are not allowed. Finally the solution of the 
model is discussed with the help of a numerical example. 

Key words & Phrases : Demand, EOQ, Inventory, finite replenishment, L 2 -system. 

AMS Classification No.: 90-XX. 


INTRODUCTION 

Most of the inventory models are developed for a single warehouse. Due 
to large stock and limited capacity of existing storage facility OW, an additional 
storage facility RW which is located away from OW is hired to store excess 
items. In practice large stock attracts the management to either an attractive 
price discount for bulk purchase or the management resorts this practice where 
the acquisition costs are higher than the holding cost in R W or due to existing 
storage capacity limitations. 

One warehouse, inventory storage facility is known as Z,- system whereas 
a two warehouse inventory storage facility is known as Z 2 - system. A good number 
of authors have considered Z 2 - system under different assumptions. Hartley [5] 
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discussed this type of model neglecting the cost of transporting a unit from RW 
to OW. Sarma [11] extended it by considering the transportation cost of a unit 
from RW to OW with AT-release rule and improved the working of the model [12, 
13]. Dave [2] also improved the working of the model given by Sarma [11]. 
Murdeshwar and Sathe [7] extended the model given by Sarma [11] by 
considering finite production rate. Gupta and Vrat [4] developed a single level 
inventory model with stock dependent consumption rate. Mendus and Antony 
[8] developed a two warehouse inventory model by considering stock dependent 
consumption rate and infinite replenishment rate. 

In this paper a L 2 - system with stock dependent consumption rate with 
finite rate of replenishment with K- release rule (KRR) is considered, as large 
stock attracts customers particularly in consumer goods. Realistic assumptions 
given in this paper make the model more practical. 

ASSUMPTIONS AND NOTATIONS 

The following assumptions and notations are used for the development of 
the model. 

(i) D is the demand rate and is taken as D = a + \Q*(a,\,8 are constant). 

(ii) Q is the highest inventory level. 

(iii) The items of RW are transferred to OW in V shipment of which K <W 
units are transported in each shipment. 

(iv) Replenishment/Production rate is finite. Let'/?' be the finite replenishment 
rate per time unit, which is known, constant and uniform during the period 
under consideration. 

(v) Storage capacity of OW is W and that of /? IT is unlimited. 

(vi) Let p = 1 - D / /? and W t = Jf7 P. 

(vii) Lead time is zero. 

(viii) Set up cost is fixed and taken as A. 

(ix) The transportation cost of K units from R W to OW is C t at a time, which is 
constant over time. 
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(x) The holding cost/unit in OW is H and that in AIFis F and F>H. 

(xi) Shortages are not allowed. 

(xii) Consumption takes place from OW only. 

(xiii) T represents the total time period (cycle length). 

MATHEMATICAL MODEL AND ANALYSIS FOR THE SYSTEM 

We assume that a company purchases/products Q units, the finite rate of 
replenishment/production is R. During production stage the demand is satisfied 
from production and the excess quantities are transferred to warehouses as in 
Fig. (i). It is economical to fill OW first and when it is filled to its maximum 
capacity then excess production/replenishment is transferred to RW. When the 
production stops, the accumulated inventory can be used to satisfy further demand 
from OW. Due to this, inventory level in OW will fall from its maximum capacity 
W. When the stock level in OW reaches to {W-K) units, at this stage 'K units 
from RW are transported to OW to meet further demand and this process is 
repeated V times untill stock in AIT is fully exhausted. The remaining (W-K) 
units in OW are used again at this stage. Also noting that after production stops, 
the inventory level in OW is 'W and that in RW is z = (Qfi-W) units. The 
inventory level is shown in figure (ii). 
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Case (i) (Before Consumption) 

Initially, we have the total average cost as 

' Q 2 2 

Where D = a + \Q S 


( 1 ) 


For optimum value Q * of Q, which minimize total average cost is obtained by 
solving 


dC(Q) I (IQ - 0 which gives Q * 



(If ?L = 0) 


As p 1 i.e. for infinity rate of replenishment 


Q* = — which is same as Classical EOQ formula. 

Case (ii) (During Consumption) 

From figure (1). It is clear that during production stage. 

Total inventory in OW = $W@I R-$W? HR 

Total holding cost of these items = H\$!>W X QI R-$W* /2/?] (2) 

Total inventory in RW = $(Q-W t ) 2 /2R ■ 

Total holding cost of these items = F^>{Q-W^ HR . 

Further as soon as the production stops [as in Fig. (ii)] 
there are W units in OW and $(Q- W { ) units in RW. 

Now using KRR, we have 


( 3 ) 
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Inventory units in RW are given by 
= t k [Z + (Z - K) + (Z - 2K)+ .+Z - (n - 1)*] 

2 

where t k = K/ D, time taken for the consumption of K units. 

Total holding cost of these items in RW = F$k(Q-W^){n + \)l2D (4) 

Cost of transporting K units from 

RWXo OW is given by =nC, =$(Q-W x )C,IK (5) 

When K units are drawn from R ITin each shipment more are carried in OW for a 
period of t k and hence holding cost of these items is given by KHtJ2. 

Since there are such shipment and taking the initial K units into account, we 
have The total holding cost of these items 

= {(n + \)HKt k } / 2 = {(« + \)HK 2 } / 2Z) (6) 

A quantity of ( W-K) units is kept unused in OW for a period of t w _ k =(n + \)t k 
and during usage the average inventory in OW is given by (W-K)/2 units and for 
a period (*-/„,_*). Hence inventory holding cost in OW of these items is given 
by 


H{K(W-K)(n + \)/D + (W-K) 2 /2D} (7) 

The fixed order cost per order is A. Hence the total cost of inventory for the 
system using (2) to (7) is given by 

C= A + H{$W { QIR-$W? I2R} + F$(Q-W X ) 2 UR + FK$(Q-W { ) («+l)/2D 
+nCt + (n + 1 )HK 2 / 2D + h[k(W-K)( n +1 )/D + (W-K) 2 /2D] 
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Average inventory cost 
C(Q,K) = C/T 

Also we have T=Q/D, Z = $(Q-W,), n = Z / K andD = a + \Q & , 

C(Q,K)= A(a + XQ 6 )/Q + pFQ/2 + pW ] 2 (F-H)/2Q-pKW l (F-H)/2Q 

+PK(F-H)/2-VW l (F-H) + pC l (a + XQ 6 )(Q-W i )/QK 

( 8 ) 

For optimum value of Q and K , substituting 6C(Q,K)/bQ = 0 and 
8C(Q,K)/8K = 0 bC(Q,K)!bQ = 0 implies 


V>+i 


^Q£{8-»;(8-i)} 


+ 2^(8-l)0 8 +pF0 2 


= 2Aa+$W?{F-H)-$KW x (F-H)-2C$WplK (9) 

which can not be solved for Q analytically and can be solved with the help of 
Newton Raphson method. 

Also if X = O 

\_ 

Q* =[2Aa/^F + W l 2 (\-H/F)-W l K(\-H/F)-2C l W i a/KFf (10) 

6C(Q,K)/bK = 0 implies 



2C,(a + X0 8 ) 
__ 


2C,a 

F-H 


(If A. * 0) 


(11) 


Also from (11) and (10) we have 
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0 * = 


2 Aa 


I P F — 


COST REDUCTION DUE TO ^-RELEASE RULE 


The unit cost of KRR is C' =C,/ K- Suppose the unit cost of transportation 

C* without bulk release rule. Bulk transportation is economical only if C* > C’. 
Hence without KRR, we have the cost function 

C{Q) = A(a+^)/Q+VFQ/2+VW?(F-rf)/2Q-W{F-tf) 


+ VC;(a + XQ*)-C;VW l (a + XQ s )/Q (13) 

The optimum value Q' of Q which minimize (13) is obtained by solving 
dC{Q)/dQ = 0 
This implies 

pc;X80 8+1 + 0 8 [/4X(8 -1) - Cfp»^(8 -1)] + p FQ 2 = 2Aa + $W?{F- H)- 2C^p^a 

Equation (14) can not be solved for Q analytic and can be solved with the 
help of Newton Raphson method for obtaining optimal value of Q. 

If X = 0 then equation (14) reduces to 

i 

Ql=[2Aa./pF + W, 2 (F-H)/F-2C;W ] a/F] 2 (15) 

KRR is economical if [ C(Q ) - C(Q, K)] > 0 

This implies (C, -C',)>(F-H)K/2a (16) 

must be satisfied. 

Thus for a given situation inequality (16) is satisfied than KRR will be 
economical. 
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Relation between EOQ with and without KRR. 






( 17 ) 


In general Q * exceeds Q. 

SPECIAL CASES 

Case (i) If (3 -> 1 i.e. rate of replenishment is infinite. 

Then equation (10) for optimum value of Q. becomes 

\ 

Q*=[2Aa/F+W 2 (\-H/F)-WK(\-H/F)-2C'lVoi/KFf (18) 

Equation (18) is the same as equation (9) as that given by J. Mendus and Antony 

( 8 ). 

Case (ii) If Demand D = a and p —>■ 1. 

Then equations (10), (11) and (12) of our Model reduces respectively to 
the equations (6) and (7) of Model (11). 

Case (iii) When F is large, Q tends to W v 

Further when C, -> 0 Then Q * = [2Aa /^F + W 2 (l-H/ F)p 

In this case bulk transportation is not necessary. 

NUMERICAL EXAMPLE 

For a given inventory system, suppose we have H = 2, F = 4, A = 100, W = 100, 
C, =2, C,‘ = 0.15, a = 2000, \ = 0,R = 6000, D = 2000Thenp = 1 -D/= 2/3, 
W, = JF/p = 150. 

(a) With KRR, we have K ’ = 63, Q' = 389 and C(Q, K ) = Rs. 923.17. 

(b) Without KRR. Q = 372; C(Q) = Rs. 993.31 
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Also when k = 63, Q = C, IK = 0.0137 and C= 0.15 clearly C,* > C\. 
Also it is clear Q* > Q'. 


CONCLUDING REMARKS 

In the present paper, we have formulated and solved a deterministic 
inventory model for two warehouses considering finite rate of replenishment 
and stock dependent consumption rate. Profitability of KRR has been desired. 
Also optimum batch size 'AT is independent of production rate and capacity of 
own warehouse, but still from the example it is clear that KRR reduces cost 
considerably. Further it is worth remembering that one should consider the 
possibility of RW whenever lot size indicated by the classical EOQ model is 
greater than the capacity of own warehouse. This model further can be extended 
by considering shortages, as shortages is an unavoidable characteristic of most 
of the inventory systems implemented in practice. 
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OSCILLATION THEOREMS FOR FIRST ORDER 
NONLINEAR FUNCTIONAL DIFFERENTIAL EQUATIONS 

OF NEUTRAL TYPE 

H. A. SHAOUL* 

(Received 27.07.2002) 

INTRODUCTION 

In this paper we are concerned with the oscillatory behaviour of first order 
neutral functional differential equation of the form 

^(x(/) + c(0-v[x(0]) + 5^(0/(x[g(0]) = 0, (1.1; S) 

where 5 = ±l,T,c,g,g:[/ 0 oo) R = (t 0 x) and/:/?-> R are continuous. The following 
conditions are always assumed to hold: 

(i) Hms(r) = =°; 

(ii) )is strictly increasing for t > t 0 > 0 and Jjjjj T (*) = °°> 

(iii) xf(x) > 0 for * * 0; and 

(iv) q(t) >o and not identically zero on any ray of the form [£«>) for some T>t () . 

By a solution of Eq. (1.1;§) we mean a function x :[7’ x ,°o)->R such that 
.y(/) + c(/)a:[t(/)] is continuously differentiable and satisfies Eq. (1.1;§) for all large 
t > T x The solutions which vanish for all large t will be excluded from our consideration. 

A solution of Eq. (1.1 ;g) is said to be oscillatory if it has an infinite sequence of 
zeros tending to infinity, otherwise, a solution is said to be nonoscillatory. Eq.(l.l;§) 
is said to be oscillatory if all its solutions are oscillatory. 

The oscillator behavior of first order functional differential equations of neutral 
type has been intensively studied in the literature in recent years. Most of the literature, 
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however, is concerned with linear equations with constant coefficients and constant 
deviations for which oscillation criteria are given through analysis of the associated 
characteristic equations, and the references cited therein, and very little is known 
about the oscillation of Eq. (1.1 ;8) whose coefficients and deviations are variable. 

It seems that the known oscillation criteria fail to describe the oscillatory behavior 
of Eq. (1.1 ;6) when/is not a monotonic function. 

In this paper, we deal with the oscillation problem of nonlinear Eq. (1.1 ;g) with 
variable coefficients and deviations and with/not a monotonic function. We establish 
some new sufficient conditions for the oscillation of all/and or bounded solutions of 
Eq. (1.1 ;§) where the function c (t) is either positive or negative for t>t 0 . We also 
present some new sufficient conditions for the oscillation of the forced neutral equation. 

(1 -2;8) ~f t (*(0 + c(0*[x(/)]) + 5q(t)f(x[g(t)]) = <?(/), 


where e:[t 0 ccj -*R is continuous. 

Examples are inserted in the text to illustrate the relevance of our theorems. 

OSCILLATION OF EQ. (1-1; 8) 

The following notations are needed in the sequel: 

^,=K-OU[t 0 oo) if t o >0 

= (-°o,0) U (0,oo) if t o = 0 

C(R) = {/:/?-> R\f is continuous and x f(x) > 0 for * * o} • 

C b [R, ) = { feC(R): f is of bounded variation on any interval [a, 6] c R, } 

Suppose t o >0 and /eC(/?).Then ftC B [R, ) if and only if 

f(x) = G{x)H(x) for all x e R, 


where 
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G: R, -> (O,oo) is nondecreasing on (- 00 ,-f 0 ) and nonicreasing on (/ 0 00 ) and 
H : R, —>R and nondecreasing on R, . 

We assume that the functions c and x are subject to one of the following 


conditions: 

0 < c(f) < a, < 1 x(/)>t; and (1.3) 

1 <b { <c(t)<b 2 x(f)</; (1.4) 

0 < c(t) < a, < l x(r) < r; and (1.5) 

1 < b, < c(/) <, b 2 x(t)<t. (1.6) 

We let 


a = 1-a, and 


b = 


V-l 

b,b 2 


The main results of this section now follow: 


Theorem 1.1 Suppose that /fC( >0 and let the functions G and//be a pair 

of continuous components of F, //being the nondecreasing one. Then the following 
conclusions hold. 

(I) Eq. (l.l;s) is oscillatory, ifcondition (1.3) holds and for every positive constant 
c< 1, g (t) < t for t > t 0 and all large t, the first order delay equation 

/(/) + « 9 (0G(c)//(v[g(/)]) = 0 (1.7) 

is oscillatory. 

(II) Eq. (1.1; 1) is oscillatory, if condition (1.4) holds and for every positive constant 
c < 1, x _1 °g(t) <t for t > t 0 (where t -‘ denotes the inverse function of x(t)) and all 
large t, the first order delay equation 

V(t) + bq{t)G(c)H(x{x-' og(0]) = 0 (1.8) 

is oscillatory. 

Proof: Let x(t) be a nonoscillatory solution of Eq. (1.1; 1). We may assume that x(t) 
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is eventually positive, since a parallel argument holds if jc(/) is eventually negative. 
Define 

y(/) = *(/) + c(0*[ x (0]> (1.9) 

then y(r) >0 for t > t 0 and Eq. (1.1 ;5) implies that y(t) is a decreasing function for 
all large t. Using this fact, we see from (1.9) that if condition (1.3) holds, then 

*(0=;K0-c(0+(0] 

= y(t) - c(0[r[^(0] - c[t(0]*[t ° t(/)]] 

>{\-c(t))y(t) 

> (l-a,)y(t), 
or 

x(t) > ay(t), for all large t, (110) 


and if condition (1.4) holds, then 



i i 

'y 

V 

'oT-'(0] 


[V‘(oj 

i 

H 

i 

K C 

y 

ox-(0] 



cTt ’«T-’(/)l-l , , 




V b \ b 2 J 


K X "'W]. 


or 
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x(t)>b >-[x _l (/)], for all large/. (1.11) 

(A) Suppose that condition (1.3) holds. Then from (1.10) and the fact that >>'(/) s 0 
for all large t, it follows that there exist a T>t 0 and a positive constant c* < 1 such 
that for all t > j 

( 1 . 12 ) 

and 

*k(O]^0>{g(O]- (1-13) 

Using (1.12) and (1.13) in Eq. (1.1:5), we have 

o=y(o+^(0/W^)]) 

= y(r) + ^(/)G(x[g(0])//(x[g(/)]) 

> y (r) + G(c*)9(0//(ay[g(/)]), 

or 

vt/(/) + aG(c*)<7(/)//(w[^f)])<0 for t> T, (1-14) 

where w(t) = ay(t). It is easy to check that for all t>T 

oo 

w(/) > aG(c)j q(s)H(\\{g{s)])ds. 

I 

The function w (t) is positive and strictly decreasing for t>T. 

(B) Suppose that condition (1.4) holds. Then from (1.11) and the fact that y (t) < 0 

for all large t, there exist a positive constant c * < \ and aT £ t 0 such that (1.12) holds 
and 

x[g{t)]tb y[x"'og(0] for t >f. (1.15) 

Now, using (1.12) and (1.15) in Eq. (1.1; 6), we have, 

/(0 + G(c‘)q{t)H(by[x-' o g(t)]) S 0 for 1 2> f , 
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or 

u' , (/) + />G(c')^(/)//(w[T , og(/)])^0 for/^r*. (1.16) 

Where w ( t ) = by (/)• 

Proceeding as in the above case, we obtain the desired contradiction. 

This completes the proof. 

For any continuous function, a:[/ 0 oo) -» R with jin) a(r) = oo, 
we define 

R(a) = {r e[r 0 ,oo]:r 0 < a(r) < 1 }. 

Also, we assume that 

//(*).*£« * 1*1“ for x >0 and cx>0. (1.17) 

As an application of Theorem 1.1, we get the following oscillation criterion for 

Eq.(U;5). 

Corollary 1.1 Suppose that/eC(/f, u ),f 0 s: 0 and let the functions G and H form a 
pair of continuous components of F,H being the nondecreasing one. Then, 

(I) Eq. (1.1 ;5) is oscillatory if conditions (1.3) and (1.17) hold and one of the following 
conditions (^,) or (5,) holds: 

(Al) a = 1, and 

g(/)<r and g(t) is nondecreasing for t>t 0 \ (1.18) 

and for every positive constant c <1, 

im inf } q{s)ds > • (119) 

l) aeG ( c ) 


(B,) a < 1, and 

jq(s)ds = «. (1-20) 

«f) 


(II) Eq. (1.1 ;g) is oscillatory if conditions (1.4) and (1.17) hold and one of the following 
conditions (C,) or (D,) holds: 
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(C,) a = 1, and 

a(0 = t' 1 °g(t) z t anda(t) is nondecreasing for t>/ 0 and (1.21) 

for every positive constant c<l, 


t 

lim inf J <l( s ) ds 

<-+<*> * x 
o(t) 


beG(c). 


( 1 . 22 ) 


(».) 


a <1 and 


| q(s)ds = oo. 
*(' '*«) ■ 


(1.23) 


Proof: Let x (t) be an eventually positive solution of Eq. (1.1 ;5). As in the proof of 
Theorem 1.1, we obtain the inequalities (1.14) and (1.17). Now, we consider the 
inequality (1.14). 

(A,) From (1.17), inequality (1.14) takes the form 


W (/) + aG(c')q(t)w[g(tj] < 0 for t > T x > T. (1.24) 

In view of condition (1.19), it follows that inequality (1.24) cannot have an eventually 
positive solution. This contradicts the fact that w{t)-ay{t)> 0 eventually. 


(B,) Since w[g(/)] > w(/) for all r e R (g)n[7\a>), we have 


J 

*(g)n[7>] 


q(s)ds < 


1 

aG(c*) 


»(r) 

| v" a (/v 


0 


aG(c’) 


w'-°(r) 

1 - a 


< 


which contradicts (1.20) 

Next, we consider the inequality (1.16). By a similar argument, as above, we obtain 
the desired conclusion. 

The following example is illustrative. 

Example 1.1 Each of the following neutral equations 
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d_( 

dt 


x{t) + -x[t + 2Ti\\ + -e 


i.2 t 


2[l-3«] 


x[r-3n] = 0,t > 0, 


and 


-(x{t) + 2x[t - 2n]) + 3e |sin,| e~^ , ~ 5 ’ l ®.x[r - 5n] = 0,/ > 0, 


(1.25) 


(1.26) 


has an oscillatory solution x(t) =sint. 

All conditions of corollary 1.1 (I-A,) and (II-C,) are satisfied for Eq. (1.25) and Eq. 
(1.26) respectively. 

BOUNDED OSCILLATION OFEQ.(l.l; 8) 

In this section we are concerned with the oscillatory behavior of all bounded 
solutions of Eq. (1.1; 8) when condition (1.5) or (1.6) holds. Also we study the 
oscillatory behavior of Eq. (1.1) where c ( t ) is negative. 

Theorem 1.2 Suppose that/ e C (/?, o ),^ 0 > 0 and let the functions G and H form a 
pair of continuous components of F, //being the nondecreasing one. Then, 

(I) All bounded solutions of Eq. (1.1;-1) are oscillatory if condition (1.5) holds 
and 

one of the following conditions (A) or (B) is satisfied: 

(A) 

\q{s)ds = oo ( 1 . 27 ) 

0 

(B) Condition (1.16) holds with a = 1, 


g (t) > t and g (t) is nondecreasing for t > / 0 , and for every 

positive constant c > l, 

*(') j 

liminf |cr(jWj>-7—. 

'->• \ W aeG(c) 


(1.28) 

(1.29) 


(II) All bounded solutions of Eq. (1.1; 8) are oscillatory if condition (1.6) holds 
and one of the following conditions (C) or (D) is satisfied 

(C) Condition (1.27) holds; 
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(D) Condition (1.17) holds with a= 1, 


a(t) = x 1 o g[t) > t and cr(f) is nondecreasing for t > t 0 , and for every (1.30) 
positive constant c > 1, 


a(i) 

liminf J^(a)(/.v> 


_ 1 

beG(c)' 


(1.31) 


Proof: Let x ( t) be a bounded and eventually positive solution of Eq. (1.1; 8 ).Define 
y(0 by (1.9). Theny (/) is eventually positive. From Eq. (1,1 ;§) it follows that y (t) is 
nondecreasing for all large t, and it follows that (1.10) or (1.11) is satisfied according 
as condition (1.5) or (1.6) holds. Since x (/) is bounded andy (/) is nondecrasing for 
/ > t 0 , there exist positive constants (3, > 1 and P 2 > 0 and a T>t 0 such that for t > j 

*[g(0]^Pi and l^kCOj > (1.32) 

and both (1.13) and (1.15) hold. 

(1) Assume that condition (1.5) holds. Inview of (1.13) and (1.32), it follows from 
Eq. (1.1-1) 

for t >T. (1.33) 


lntegrating(1.33) from 7"to l, we obtain 


oo > y(t ) > y(T) + G(p,)//(np 2 )j q(s)ds 

T 

— ) oo as t —^ oo, 

a contradiction. 

Next, using (1.17) with a=l and (1.32) in Eq. (1.1 ;8) we have 


y(0-aG(c,) 9 (/)y[g(/)]>0 for t > T. (1-34) 

But, inview of condition (1.29), it follows that inequality (1.31) cannot have an eventually 
positive solution. This contradicts the fact thaty (t) is eventually positive. 

(II) Assume that condition (1.5) holds. Inview of (1.15) and (1.32), it follows from 
Eq. (1.1 ;§) that 

/(/)^G(P,)( 7 (/)//(^[a(/)]) for I'ZT. 
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Proceeding as in the proof of the above case, we obtain the desired conclusion. This 
completes the proof. 

For illustration we consider the following example: 

Example: 1.2 Consider the neutral equations 


dj 

dt \ 


x(/) + -^x[/-c]j = (2 + e cosht' 


J£L 


= 0 (> 0 , 


and 


^(*(/) + 2x[/ + <r]) = (2f?"° +e" 2 °) cosh (e'* 2 °) 

seek (*[/ + 2a])x[/ + 2a] = 0, 
where a. is a positive constant, a > 2. 


(1.35) 


(1.36) 
t> 0 , 


All conditions of Theorem 1.2 (1-B) and(II-D) are satisfied for Eq. (1.35) and 
(1.36) respectively, and hence all bounded solutions of Eq. (1.35) and Eq. (1.36) are 
oscillatory. 

We note that each of the equation (1.35) and (1.36) has an unbounded 
nonoscillatory solution *(/) = e'. 

Next, we present two criteria for the oscillation of all bounded solutions of 
Eq.(l. 1 ;5 ),§ = ±1, when the function c (t) is assumed to be negative. 

We assume that 

0 < -c(t) = p(/) < P 0 for some positive constant P 0 . (137) 

Theorems. Let/e C (/?, o ), t 0 > 0 and conditions (1.18) and (1.37) hold and suppose 

that the functions G and H form a pair of continuous components of F,H being the 
nondecreasing one. Suppose that for every positive constant c and all large t the delay 
first orders equation 


z , {t) + q(t)G(c)H(z[g{t)]) = 0 


(1.38) 
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is oscillatory. Then every bounded solution of Eq. (1:1; 1) is oscillatory provided that 
either condition (1.27) is satisfied or in addition to conditions (1.17) with a = 1 and 
(1.30) the condition 


a(r) 


liminf | q(s)ds : 


1 


eMC) 


(1.39) 


is satisfied. 

Proof: Let .*(/) be abounded and eventually positive solution of Eq. (1.1;1). 

Define 

>{<) = *(<)-«<)*[*«]■ (1.40) 

Then, inview of Eq. (1.1:5). 

y'(t) = -q (/)/(*[g(f)]) 5 0 eventually and hence y (t) is nonicreasing for all 
large, t Thus,y(t) is of one sign for all large t, say t> r, > / 0 . 

We consider the following two cases: 

Case 1 : Lety ( t ) > 0 for t > /, . There exist a positive constant y and a t 2 >t l such 
that for allt>/ 2 . 

,[g(0]< T (1.41) 


and 


(1-42) *[g(0]^>fe(0]. 


Consequently Eq. (1.1; 1) reduces to 


/(0+ 50 for / £ t 2 . 

Proceeding as in the proof of Theorem 1.1, we obtain the desired contradiction. 

Case 2: Letw(0<0forr>r,. Using (1.41), we get 
0 <-v(0 = u (r) =p(/)x[t(r)]-x(r) 
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< P(/)+(/)] 

5 p o (0*[t( 0] for t > t 2 

and here, 

*(')-r u [ T ~'(')] for ******* 

Po 

From Eq. (1.1 ;5) it follows that 

w '(0 * ^-G(Y) 9 (/)//(w[a-' o g(/)]) for t > t 
Po 

where u (t)=(3 0 w(f). 

The rest of the proof is similar to that of Theorem 1.2 and hence is omitted. This 
completes the proof. 

Theorem 1.4 Let/e t 0 >0, conditions (1.21) and (1.37) hold, and let the 

functions G and H form a pair of continuous components of F,H being the nondecreasing 
one. If for every positive constant c and all large t the delay first order equation. 

z ' ( / ) + fT G ( C )^( / ) // ( z [ T °s( / )]) = 0 (1.43) 

Po 

is oscillatory, then every bounded solution of Eq. (1.1 ;-l) is oscillatory provided that 
either condition (1.27) holds, or in addition to condition (1.17) with a = 1 and (1.28) 
the condition. 


*(») 

liminf |g(5)rf5: 


1 

~7T\ holds. 
eG(c) 


(1.44) 


Proof: Let x (/) be a bounded an eventually positive solution of Eq. (1.1; 1). Define 
y(t) by (1.40). Then, it follows thaty (/) is nondecreasing and is of one sign for all 
large t. Now, we consider the following two cases: 


Case 1: Suppose thaty (/) >0 for t £ f, > t 0 . Then there exist a t 2 £ /, and a positive 
constant y such that 

*U(0] ^ Y and *[g(0] > w[g(/)], for t > t 2 
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and hence 


y(t)> G{ y to(/)ff(y[g(0]) for t > t 2 . 

The rest of the proof is similar to that of Theorems 1.2 (I) and hence is omitted. 
Case 2. Suppose thaty (r) < 0 for t > r, > t 0 . Set 

0 >-y (0 = u (0 

and proceeds as in the proof of Theorem 1.3-case 2 to obtain 

v ’(0 + ^- G ( y )^( 0 //( v [ x ' , o g (/)])> 0 for ,>, 2 

Po 

where u(t) = P 0 v(t). We use an argument similar to the one used in the proof of 
Theorem 1.1 (II) to complete the proof. 

The following examples are illustrative: 

Example 1.3 The first order neutral equation 


j t W/)-yx[r + 2n]) + 


i y Y !+('-*) 

t 2 (t + 2nf 


A 


t -71 


(1.45) 


1 + jc 2 [/-it] 


= 0, t>n+y. 


has abounded nonoscillatory solution x(t) = t- Here y is any positive constant. 
All conditions of Theorem l .3 are satisfied the equation 




1 1 

1 y ] 

1 + (t - n) 2 ' 

V 1 + C 2 J 

{*' (/ + 2a ) 2 j 

•-* J 


y(‘- 7t ) = °, 


for y > 0 and t > %, 

has a positive solution y (t) with y (0 = 0 for every positive constant c. 


Example 1.4 Consider the neutral equation 
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(1.46) 


where y is any positive constant. It is easy to check that for any positive constant c, 
the equation 




is oscillatory (see Corollary 1.1). Therefore, all conditions of Theorem 1.3 are satisfied 
and hence all bounded solutions of Eq. (1.46) are oscillatory. Note that Eq. (1.46) 
has an unbounded nonoscillatory solution x(t) = t. 


Example 1.5 The first order neutral equation 


7 , w ') - 2 * - 1 ])= (- V 1 ) 1+c0 ' i>0 < 1 - 47 > 

has a bounded nonoscillatory solution *(/) = e'. Here x(r) = t - l,g(r) = t-2 
and hence t _l o g(/) = t -1 > /. 

All conditions of Theorem 1.4 are satisfied, the delay equation, 

y (,)+ 2(7+7) (^7^ ](i+H< -1]=°. < > ° 

has a positive solution y(t) with y(t) = 0. 

1.4. Oscillation of Eq. (1.2; 5) 

In this section we establish some sufficient conditions under which Eq. (1.2;5), 
§ = + 1 is oscillatory. 

We list the assumptions which are needed below: 


o 

V 

II 

V 

II 

«**». 

1 

Q 

a is a positive constant; 

(1.48) 

c(t) = b x - b 2 > 1,t(/) = / + o, 

a is a positive constant. 

(1.49) 
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Also, we assume the following conditions on the forcing term: 


There exists a differentiable and oscillatory function 

ti : [/ 0 ,°°] -> R Such that (/) = e (t), (1.50) 

and 

r| (0 is periodic of period a i.e. (1.51) 

q(r±cr) = r|(f), where o is defined as in (1.49) and (1.50). 

Theorem 1.5. Let condition (1.50) hold, c ( t ) = a, a is a nonnegative constant and 
x(r) = t + 8<r, 8 = ±1, and a is a positive constant. If 

(1.52) limsupr|(f) = +oo and liminfii(/) = -oo, 

then, 

(i) All bounded solutions of Eq. (1.2; 5) are oscillatory, 

(ii) Eq. (1.2; 5) is oscillatory. 

Proof: Let *(/) be abounded an eventually positive of Eq.(1.2;$) and x(f)be an 
eventually positive solution of Eq. (1.2; 5) Define 

i y(/) = jc(/) + ax[/+8 l CT] and x(t) = y(t)-r\(t), (1.53) 

where 8, = ±1. 

Then Eq. (1.3;) take the form 

z' (t) = -8g(0/(jr[g(/)]). 0-54) 

Next, we consider the following two cases: 

Case 1.8, = -1. It follows that z ( t ) is an increasing function, say on [*,») for some 
/, ^ t 0 > 0. We show that z(t) > 0 for t> T > /,. 

If not, then z (/) < 0 for / £ T > t L Hence 

y(0-ri(/) > 0 i.e 0 >y(t) < n(0 for fz J, 
which is a contradiction, since t) (/) is oscillatory. Thus, we have 
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z(t) > 0 and z'(0)>0 for t>T. 


(1.55) 


Taking lim superior iny (/) >t| (r) as t -» oo,we obtain 
limsup>>(/) > limsupr|(r) = oo 

which contradicts the fact thaty (/) is bounded. 

Case 2.8, =1. From (1.54), it follows that z (t) is an eventually of one sign. As in the 
proof of case 1, above, we see that 

z(t ) >0 and z'(t)< 0 for t> T. (1 -56) 

Since z (0 + r|(r) = y(t) > 0, we have 
z(t) > -r\(t), from which it follows that 
lim sup z(t) £ limsup(-r|(f)) = -lim ifn r|(/) = oo 

V ' t-4-r. ' v /f !-+T) v 7 

which contradicts the fact that z (t) is bounded from above.This completes the proof. 

For illustration we provide the following two examples: 

Example 1.6 Consider the forced neutral equation 



+ ax[t + 8a]) - 


(l + ae‘ 8o )e-/ 

(l-e-^y 




1+JC 


m u , 


= t.cost+sin t, t>0, 


(1.57) 


where 8 = ±1, a and a are nonnegative real numbers, a is a positive constant, g: 
[/ 0 ooj — >/? is continuous and Jimg(/) = oo.if we choose r\(t) = t.sint, then all the 
hypotheses of Theorem 1.5 (i) are satisfied and hence every bounded solution of Eq. 
(1.57) is oscillatoiy, It is easy to veri fy that the corresponding unforced neutral equation, 
namely 


—(x{t) + ax | 


t'+H)- 


( l+ 


ae 






(1.58) 
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j^yg^[g{t)]=o,t>o. 


has a bounded nonoscillatory solution x(t) = l - e ~'. 


Example 1.7 Consider the forced neutral equation 
x(t) + ax\t + 8o]) + (l + ae ria ) 


(159) 


.(H-e*') 




a ^ 


1 + 4*0 


s £ w *[g0)] = e‘(sin t+ cost), 


where 8 = ±1, a and a are nonnegative constants, a > 0, and 
g: (/„«) R is continuous and = °°Here, we choose q(r) = e‘ sin t and 

find that all the conditions of Theorem 1.5 (ii) are fulfiled. Thus, Eq. (1.59) is oscillatory 
we also note that the corresponding unforced neutral equation 


dt 


(x + ax[t + 8cr]) + (l + ae ' 8o )(l + e g{,) ) 


(1.60) 


1 ' 




sgn x[g(/)] = 0, t >0, 


has a nonoscillatory solution x(f) = e". 


Theorem 1.6 Let/ e c (/? (| ),/ 0 ^0 and conditions (1.50) and (1.51) hold. 

Suppose that G and H form a pair of continuous components of F, H being the 
nondecreasing one. Then the following conclusions hold: 

(I) All bounded solutions of Eq. (1.2; 5) are oscillatory if condition (1.48) holds and 
either 

(A) Condition (1.27) holds; or 

(B) In addition to condition (1.7) with « = 1 and (1.28) the condition (1.29) holds 
for every constant c > 1 
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(II) All bounded solutions of Eq. (1.2; 5) are oscillatory if condition (1.49) holds and 
either: 

(C) Condition (1.27) holds: or 

(D) In addition to conditions (1.17) with a = 1 and (1 -30) the condition (1.31) 
holds for every constant c > 1. 

Proof: Let *(/) be a bounded and eventually positive solution of Eq. (1.2; 5). As in 
the proof of Theorem 1.5, we define the functions jy (t) and z (t) as in (1.53) and 
obtain 

z{t) > 0 and f (t) > 0 eventually. (1.61) 

Next, we consider: 

(I) Let condition (1.48) hold. From (1.50), (1.51), (1.53) with 8, = -1 and 
(1.61), we see that 

x(f) = z(t)+ r|(t) - a l x[t - o] 

= z(t) + r|(t) - n, [z(t - a) + r|[t - a] - a i x[t - 2a]] 

> a(z(t) + q(/)) for all large t, 

where a=l-a { . Next by (1.51), there exists a T> t 0 (say) such that for all / > T 
x{t)>a(z{t) + r\{T))=w{t). 

Clearly az' (t) = W (/) and w (t) > 0 for all t > T. Since x(t) is bounded, there exist 
a T X >T and a positive constant y such that for / > T x 

x[g[/]]<y and x[g(t)] > w[g(/)]. (1.62) 

Thus, W (t) > aq(t)G{y )H( w[g{/)]) for all t>T r 

The rest of the proof is similar to the proof of Theorem 1.2 (I) and hence is omitted. 

(II) Let condition (1.49) hold. From (1.50), (1.53) with 8, = 1 and (1.61), we obtain 

x(t) = Uz[t - a] + r)[t - o]-x[t - a]) 
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= -J- (z{t - a] + rj[/ ~ a]) - (z[t - 2 a] + v[t - 2a] - x[t - 2 a]) 

0 \ 0J 

v °\ J 

= b(z[t - a] + ri[/ - a]) for all large t, 

i _ — 1 

where ° ~ ,2 • 
b \ 

Once again by (1.51), there exists a T * >t ] + a for some 
r, > t 0 Such that for / > j* we have 

x{t) * b (z[t - a] + Ti[r* - a]) = v[/-a], 

and bz' (/) = V (/) and v(r) > 0 for t>T* ■ Since x(t) is bounded there exist a constant 
y, > 1 and a T' > T* such that 

*[$(')] 5 Y'i and x[g(/)] > v[g(/)-a] for all t > T* (1.63) 

Thus V ( t) > bG( y, )q(t)H(v[g{t) - a]) for t > T *. 

Proceeding as in the proof of Theorem 1.2 (II), we obtain the desired conclusion and 
the proof is now complete. 

For illustration we consider the following two examples: 

Example 1.8 Consider the forced neutral equation 

^(x(t) + cx[t - 2it5, ]) = ce-'V 1 ' 1 *"') (1,64;5,) 

(K«mD sgn x[g(t)]+co5r, t >0 
where a > 0,5, = ±1 and c is a positive constant such that 

(i) 8,=1 if 0<c<l;and 

(ii) 8,=-l if Ol, 

g:[t Q 00 ) R is continuous and g (r) = 00 . 
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Here, we choose the functions, g, H and r\ as follows 

G(x) = e' M ,H(x) = |*|“ sgn x and r)(r) = Sin t. 

Then r|(/) is oscillaory and periodic of period 2 tt. Also, 

if (/) = cost= e(t). 

All conditions of Theorem 1.6 (I-A) and (II-C) are satisfied for (i) and (ii) respectively 
and hence all bounded solutions of Eq. (1.64;8,) are oscillatory. 

i \ 3rc 

We observe that if a = l and = —, Eq. (1.64; 8,) has a bounded and 

oscillatory solution *(/) =sint. 


Example. 1.9 Consider the forced neutral equation 


—(*(/) + «[/-2*8,]) = 


f 1 + c + ccost ^ 
tanh(/ + sinf) 


(1.65;5) 


5 ec/i(^[g(/)])sinh(x:[g(/)]) + cosi, t >0 

where 8, = ±1, c and g are defined as in Eq. (1.164;8,) and (i) and (ii) in Example 1.8 
hold. Here, we let 

G(x) = sech x, H (x) = sink x and t](r) = sint. 

All conditions of Theorem 1.6 (I. A.) and (II.C) are satisfied for (i) and (ii) hold 
respectively, and hence every bounded solution of Eq. (1.65;s) is oscillatory, It is 
easy to check that ifg (t) = t, Eq.( 1.6;8,) has an unbounded nonoscillatory solution x 
x{t) = / + sint, t > 0. 


SOME GENERAL REMARKS 


1. We note that if/(x) in Eq. (l.;8,) and (1.2;8,), 8, =±1, is a nondecreasing 
function for** 0(i.e.,the components G (*) of F(x) is identically 1), then the 
boundedness requirement on the solution of Eq. (1.1.8) and Eq. (1.2;5) in 
Theorems 1.2,1.3,1.4 and 1.6 can be removed. 

2. The function f(x) in Eq. (1.2;8), 8,= ± 1 may not satisfy the linearization 
conditions. For example, we see that if 


Ax) i*r . i 

* l + x 2 


as *—>0* 


a = 1 
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-> 1 as x -> ±oo, a = 3 

—> oo as x —> ±oo, a = 3- 

The lineariztion conditions are satisfied while for the case when 

= e' w |jc| 0 ‘' -> 0 as x -> 0 and x -> ±oo, a > 1. 
x 

The linearization conditions fail. Thus, the results of this paper are applicable to a 
larger class of equations of type (1.1; 8,) and (1.2; 5,), 8, =+l. 

3. The forcing term e (t) in Eq. (1.1; 5,) and the function r| defined by (1.50) and 
(1.51) are assumed to be oscillatory and periodic and hence bounded for t > t 0 . 
We note that: 

(i) The forcing term e considered in Eq. (1.2;8,) is not necessarily small. 

(ii) From Examples 1.6 and 1.7 it is evident that the presence of a forcing term 
can generate oscillations in an otherwise nonoscillatory equation. 

Finally, it should be interesting to obtain results similar to Theorems 1.2-1.6 for 
Eq.( 1.11) and Eq. (1.2;-1) without the boundedness restriction on the solutions 
under consideration. 
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PROPAGATION OF STRONG SPHERICAL DIVERGING 
SHOCK WAVE AND CHANGE IN ENTROPY OF 
NON-UNIFORM MEDIUM 

R. P. YADAV* AND P K. GANGWAR* 

(Received 15.04.2002 and in revised form 06.09.2002) 

ABSTRACT 

Chester-Chisnell-Whitliam method has been used to study the propagation of shperical 
diverging strong shock wave in non-uniform medium. Assuming an initial density distribution 
law as p 0 = p'r"“ where p' is the density at the centre and to is a constant, analytical relation 
for shock velocity has been derived. The expressions for pressure, particle velocity, sound 
velocity, temperature and change in entropy of the medium immediately behind the shock have 
been obtained. Finally, the results accomplished here have also been compared with earlier 
results as well as with those for converging shock. It is found that change in entropy increases 
as shock propagates in the medium. This agrees with the universal law of increase of entropy. 

Key words : Diverging shock wave (Non-uniform medium), Change in entropy 

INTRODUCTION 

On account of its immense importance in Astrophysics and 
Thermodynamics, propagation of shock wave (plane, cylindrical, and spherical) 
has received considerable attention of many workers [1-10]. Considering the 
effect of overtaking disturbances, the propagation of strong diverging shock in 
non-uniform medium, the temperature variation has been computed by Yadav 
and Rana [11]. 

Yadav et a/ [12-13] have studied the propagation of spherical converging 
strong shock in uniform medium and obtained the change in entropy and 
temperature of non-uniform medium on propagation of spherical converging 
strong shock wave. 

Recently, Yadav and Gangwar [ 14] have studied the propagation of spherical 
converging strong shock ir. non-uniform medium. In the present paper CCW 
method has been used to study the propagation of spherical diverging strong 
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shock in non-uniform medium. Assuming an initial density distribution law as 
p = pV“ where p' is the density at the centre and © is a constant, analytical 
relations for the shock velocity, the pressure, the particle velocity, the sound 
velocity, the temperature and the change in entgropy of medium immediately 
behind the shock have been derived and their variations with propagation distance 
(r), constant (©) and adiabatic index (y) have been given through figures. 

Finally, the dependence of these parameters on adiabatic index (y), 
propagation distance (r) and constant (co) have been discussed. 


It is found that temperature and change in entropy of the medium increase 
as shock propagates. This agrees with earlier results [8, 13] as well as universal 

4y 


© >■ 


law of increase of entropy. For |2 + ^2y(y -l)|» can 8 e > n entropy 

increase with © and decrease with adiabatic index y and the change in entropy 

4y 


© < 


is reverse 


for {2 +1 /2v(y~l)} 


THEORY 

The equations governing the spherically symmetrical flow enclosed by the 
shock front are : 

du du 1 dp 

dt dr p dr ' ' 


(d d \ (du -2u\ 

-h U p + p H 

[dt dr) \dr r ) 


= 0 


( 2 ) 



(3) 


where, u(r,t), p(r,t) and p (r,t) denote respectively the particle velocity, the 
pressure and density at a distance r from the origin at time t, and y is the adiabatic 
index of gas. 
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Let p 0 and p 0 denote the undisturbed value of pressure and density in front 
of shock and />,, p, and m, be the values of respective quantities at any point 
immediately after passage of shock, the Rankine-Hugoniot conditions are : 


2 A/ 2 y-1 

Y+l y(y + l) 


(4) 


(Y + 1 )M 2 

P> ~ P ° (Yi) A/ 2 + 2 


= 


2a n 


Y + l 


M- — 
M 


\ 


a ~a 0 


,J[2yM 2 -(y-l)J(y - l) w! +2 ] 
(y + 1)M 


Where M = U /a Q is Mach number, U is the shock velocity, a and a 0 are 
the sound velocities in disturbed and undisturbed medium respectively. 

The temperature and change in entropy are given in MaGraw-Hill 
Encyclopedia of science and technology [3] 


T» (y + l) ! M ! 


= - (t -1)][2-I-(y-l)Af 2 ] r 

* Y-1 n (y + l) r+1 M 2t 


(5b) 


Where T and T 0 are the temperature of perturbed and unperturbed media 
respectively, AS is the change in entropy between perturbed and unperturbed 
media and R is the universal gas constant. 
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For strong shock i.e. M» 1, conditions (4) and (5) reduce to 


2p 0 t/ 2 y + 1 2 U V 2 nY-l) rr 

p = p = - L --p 0 , u = —fl = --— U 

y-1 y-1 y + 1 y+1 


r = 2y(y-l) r £/ 
^0 (y + 1) 


— 1 Inf 2y ( y ~ T ) Y If— 

y-1 n { (y + l) T+1 jUo 


For spherical diverging shocks, the characteristic form of the system of 
equations (l)-(3), i.e. the form in which each equation contains derivatives in 
only one direction in (r, t) plane, is 

, , 2pa 2 w dr 

dp + padu + -r— = 0 ( 8 ) 

[u + a) r v ’ 

It is assumed that the medium in which shock propagates is non-uniform, 
whose initial density distribution is given by 

Po=P , r~ a ( 9 ) 

where p' is the density at the centre and o is a constant. 

Substituting conditions (6) into equation (8), we get 


d U 1 4y dr 

—— +- } . —co +- } - — = 0 

U 2 + j2y/(y-l) 2 + J2y(y-l) r 


On integrating equation (10), expression for shock velocity can be written as 


U = Kr 


W2y/(y-1) 2 +V 2 y(y-l); 
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where K is the constant of integration. 

The expression for non-dimensional shock velocity, there after called as shock 
strength can be written as 


u 


= Ay- V2 r l 


1 

„ 4y 1 

W2 Y /(y-1) 

| 2+V2y(Y-l) J 


( 12 ) 


where A = ky[pj~p 0 

From equation (6), (7) and (12), the pressure, the particle velocity, the sound 
velocity, the temperature and the change in entropy are given by expressions (13)- 
(17) respectively. 


P_ 

Po 


2_^_ 
(y + 1) 


2 

m _, 4 .Y. . 1 

2+j2yl{y-l)\ 

2+V2 Y (y-I) J 



1 


1* 

11 

to 


(0 4y | 


« 0 (y+i)Vy 


a 

“o 


L 

T 

i o 


(t + 0 ' 


2^(T-l) I 
(v + 1) 2 


1 

4y ; 

_2+V2y/(y-I) 

2+V 2 r(Y-i)| 

2 

4y 

2 + W-i) 

° 2+V2y(y-0 J 


AS 

R 


1 

(y-1) 


In 


2^ 2 (y~1) T 

(r+r 


2 

L 4 Y 11 

_2+V2y/(r-1) < 

1 W 2 Y(r-oj. 


(13) 


(14) 


(15) 


(16) 


(17) 



114 PROPAGATION OF STRONG SPHERICAL DIVERGING SHOCK WAVE AND CHANGE IN 



17.S-- 
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r 


Figure I Variation of shock vriocirv with propagation distance <r) for 

adiabatic index r"U3 



Figure 2 Variation of pm sure (p/p«) with propagation distance (r) for adiabatic 
index y"1.33 
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Figure 3 Venation of particle velocity (u/ao) with propagation distance (r) for 
adiabatic index r 1.33 



-wi.y 


Figure 4 Variation of sound velocity (a/e*) with propagation distance (r) for 
adiabatic index y-1-33 
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adhtatie index r"1.33 
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Y 

Fig 7 Variation of shock velocity (L7an) with adiabatic index y at propagation 
distance r » 10 
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Fig Q Variation of pressure (p/po) with adiabatic index y at propagation 
distance r * 10 
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Fig io Variation of sound velocity (a/a«) with adiabatic index 7 at propagation 
distance r* 10 
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Fig ilVuiktioa of change in entropy (AS/R) with adiabatic index y at propagation 
distance r * 10 


rl tn 
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RESULTS & DISCUSSIONS 

Expression (12) representing, the propagation of strong spherical 
diverging shock through a non-uniform medium, shows its dependence on 
adiabatic index (y), constant (to) and propagation distance (r). It may be 
noted that for a given value of y, increase or decrease in shock velocity 
dependence on to. If to >4y /12 + ^2y(y - l)j, shock velocity increases and 
If to <4y/|2 + ^2y(y-l)} , it decreases as shock progresses. For 
to = 4y /12 + ^2y(y-l)| shock moves with constant strength. 

Equation (13)-( 15) represent the non-dimensional pressure, particle 
velocity and sound velocity immediately behind the shock of the medium as 
strong spherical diverging shock propagates, respectively. 

Taking the initial shock strength U/a 0 = 20 at r = 10, for y- 1.33 
and co = 2, the variation of pressure, particle velocity, sound velocity, 
temperature and change in entropy with propagation distance ( r ), constant 
(to) and adiabatic index (y) are numerically calculated and displayed 
through figures (1-12). 

It is found that all these parameters increase as shock propagates in 
the medium and also increase with increase in © [fig. (1-6)]. Non 
dimensional shock velocity, pressure and particle velocity decrease as 
increase in adiabatic index (y) [cf. Fig (7-9)] whereas sound velocity and 
temperature of the medium increase with increase in y . [cf. Fig (10-11)]. 

The change in entropy with propagation distance (/*) for co = 1.9, 2.0 and 
2.1 has been shown in fig. (6). It is seen that change in entropy of medium 
increases as shock propagates. This agrees with earlier results [8, 13] as 
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well as universal law of increase of entropy. The change in entropy increases 
with a) [cf. Fig. (6)] and decrease with y [cf. Fig. (12)]. It is important to 
note that for © = 4y/j2 + ^2y(y - l)j, change in entropy as well as non- 
dimensional temperature is constant and this constant depends upon adiabatic 
index (y) [cf. Fig. (16, 17)]. For values of © <4y/|2+^2y(y-l)|, the 
change in entropy of medium decreases as shock diverges, is an important 
and peculiar feature. 
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HOPF-BIFURCATION AND STABILITY OF 
A 2-PREDATORS, 1-PREY SYSTEM 

IDRISH H. ELMABRUK* 

(Received 02.12.2002) 

ABSTRACT 

This paper considers a system modelling the interaction of two competing predatory 
populations with density-dependent death rats, living exclusively on a common prey. The 
existence of interior positive equilibrium (coexistence) is proven and sufficient conditions for its 
asymptotic stability are given. Furthermore, a derivation of an easily verifiable conditions for 
the existence of a hopf-bifurcating one-parameter family periodic solutions emanating from the 
interior positive equilibrium of a general of species predator-prey system are given. The stability 
of the bifurcating periodic 3-solutions is shown, through an example, to be an open question. 

Key words: 3-species predator-prey system, hopf-bifurcation, equilibrium, periodic 
solution, stability. 

Mathematics Subject Classification (2000) : 92 BOS 

INTRODUCTION 

For a predator-prey system (see [ 1 ]) it is well known that the existence and the 
stability of a positive equilibrium is essential to investigate the existence of non-trivial 
periodic solutions. This question becomes even more interesting and important when 
two or more predatory populations compete for a single prey population. In this work, 
we have considered a system of three O.D.E.,s modelling the interaction of two 
predatory populations living exclusively on a common prey population. For both 
predatory populations, we have incorporated density-dependent death rates as 
opposed to constant death rates as considered in [3] and others, a situation that does 
not appear to be realistic. For a system with constant death rates, there may be no 
positive equilibrium existing and hence constant coexistence of the species is excluded. 
However, it is posibble that these competing predators together with their sole food 
source may coexist in a periodic fashion. In fact this does happen and the three species 
coexist in the form of state peiodicity see [ 10,11 ]. 
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In the case, however of a system with density-dependent death rates, we have 
proved the existence of the interior positive equilibrium, its asymptotic stability, and 
existence of bifurcating periodic solutions emanating from this equilibrium point. 

The organization of this paper is as follows. 

In section 2, we give the model and discuss the asymptotic stability of equilibria 
in the co-ordinate planes. 

Section 3 deals with the existence and asymptotic stability of the interior positive 
equilibrium. In section 4, we derive conditions for the existence of Hopf-bifurcating 
periodic orbits emanating from the interior positive equilibrium of a general 3-species 
system. We also show by an example that the question of stability of the bifurcating 
periodic orbit is inconclusive. 

THE MODEL 


Hypothesis on the Model >' 

A system modelling the interspecific interactions of two competing predator 
populations and living exclusively on a common prey is of the form: 

x = xf(x)-yp{x)-zq{x) 

> , = t[-^W + c,pW] (2.1.1) 

A 

z = z[-h(z) + d l q(x)], 

*(0) = x 0 > 0, y(0) = y 0 > 0, z(0) = z 0 > 0. 

where 


d_ 

dt 


> c i» d\ are positive constants known as food (biomass) conversion rates. 


We make the following assumptions, which are consistent with models of 
predator-prey system. 

(HI): j{x)\ specific growth rate of prey x,f[Q) > 0, f(x) £ 0 for x £ 0, there is 
k > 0 (the environmental carrying capacity) such that /(*)>0 on 

0 S x < k, f(k) = 0, and f(x) < 0 on x > k. 
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(H2): /?(*); the functional response of predator y with respect to prey x, and 
/>(0) = 0, p'(x) > 0 for x > 0. 

(H3): q(x) \ the functional response of predator z with respect to prey x, and 
^(0) = 0, ( 7 '(.v) > 0 for x > 0. 

(H4): q(x) = ap(x), a is a positive constant. 

(H5): g(y); density-dependent death rate of predatory, g(0)>0, and 
g'(y)>0 fory >0. 

(H6): h(z)\ density-dependent death rate of predator z, h{ 0) > 0, and h'(z)> 0 for 
z > 0. 

Asymptotic Stability of Equilibria in Co-ordinate Planes 

Clearly the system (2.1.1.) has an equilibrium at £ 0 (0,0,0) which is hyperbolic 
and unstable along the x-axis. 

By the assumption (HI), £, (fc,0,0) [a predator-fee state] is also an equilibrium 
point. The eigenvalues of the variational matrix F(£,) of the system (2.1.1.) about 
£,(&,0,0) are: 

X = kf'{k)< 0, X 2 =-g(0) + c lP {k), X 3 = -A(0) + ^#) (2.2.1) 

Remark 2.2.1. From (2.2.1.) it follows that £, (k, 0,0) is asymptotically stable along 
the x-axis. This implies that the prey population remains in a neighboured of k. For 

existence of E(x, y,0) and E(x, 0, z) as equilibria, it is necessary as well as sufficient 
that: 

-g(0) + c 1 /?(Ar) > 0 and -h(0)+d i q(k)>0 (2.2.2.) 

as it implies increase of predator populations y(/) and z(t) respectively, see [4]. Thus 
we assume that each of the predators y and z can survive on the prey x in the absence 

of its rival, that is, there exist equilibria E(x t y, 0) and £(x,0,z) satisfying: 
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xf{x)-yf{x) = o 
~g(y)+c x p{x) = o 

and 


(2.2.3.) 


xf(x)-zq{x) = 0 
-h(z) + d t q(x) = 0 

respectively. 


(2.2.4.) 


Lemma 2.2.1. If — 

A*),,.. 

E(x,y, 0) is asymptotically stable 


< 0 and h( 0) - d x q(x) > 0, then the quilibrium 


Proof: The entries of the variational matrix v{E) of the system (2.1.1.) about 

E(x,y, 0) (where simplification is done using (2.2.3.)) are: 


( -\d xf(x) 

a. i = -TT 

dx{ p(x) 

a 2l =c,yp(x)>0, 

a = 0, 


< 0, a l2 = -p(x) < 0, a l} = -q(x) < 0 

a 22 = -yg'{y) < a 23 = 0 

a J2 = 0, a 33 = ~(h( 0) -d x q(x)) < 0 


(2.2.5.) 


The characteristic equation of V[jE} is 

A, 3 + /?|A, 2 + p 2 X + p 2 — 0, (2.2.6.) 

where 


/>, =-(fl„+a 22 +fl 23 ) 

Pi = { a U a 22 +a 22 a 33 +a il a 33) - ( fl 12 a 21 + °l3 a 31 a 23 a 32 ) ' 

P 3 = ( a ll fl 23 a 32 +fl 22 a 3! a !3 + a 33 a l2°2\) ~{ a u a 22 a 2i +fl 12 fl 23 a 31 +fl 13 a 2i a 32). 


(2.2.7.) 
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From (2.2.5.) and (2.2.7.), we have 

Px > 0, p 3 > 0, and p ] p 2 - p 3 > 0. 

Thus by Routh-Hurwitz criterion, E(x,y, 0) is asymptotically stable. 


Lemma 2.2.2. If 


dx 


•*/(*) 


7v=.t 


< 0 and g(0) c,p(x) > 0, th en equilibrium 


£(x,0,z) is asymptotically stable. 


Proof : The proof is similar to that of lemma 2.2.1. 

Remark 2.2.2. There is no equilibrium in the y-z plane, for if we suppose that there 
does exist an equilibrium £*(0,y’,z’) in the y-z plane which is given by g(y*) = 0 

and h(z ) = 0, then this contradicts assumptions (H5) and (H6). Moreover, it is, indeed, 

consistent with the intuitive ecological observation that predators can not coexist while 
their sole food source is extinct. 


EXISTENCE OFINTERIOR POSITIVE EQUILIBRIUM AND ITS STABILITY 

For the system (2.2.1.) with g(y) and h(z) as contant functions, the question 
whether or not it is possible to have all three species x, y and z coexisting is a quite 
significant and interesting question both from mathematical and biological points of 
view. 


In general, no steady state x * 0, y * 0, z * 0 exists in this case. 

This is because the algebraic equations (say with Holling type II functional responses) 
xf{x)-yp{x)-zq(x) = Q 


-s + c,/>(x) = 0 
-d + d i q(x) = 0 

do not have a simultaneous solution. In the degenerate case, x, =x 2 < k, there exists 
a line segment of steady states joining £(x,y,0) and E(x,0,z). This implies, in general, 
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that the competing predators would exclude coexistence in a constant steady state 
situation. Since no constant equilibrium exists, this could occur only through a stable 
time varying solutions. In fact, this does happen and the ghree species can coexist in 
the form of a stable periodic trajectory see [ 10,11]. 

However, when both g(y) and h(z) are not constant functions, there does 
exist an interior positive equilibrium as we show now. 

Existence of Positive Equilibrium 

For the existence of positive equilibrium, we give an indirect proof. If such a 
positive equilibrium £*(x\y\z*) of the system (2.1.1.) exists then it is given by: 

x f(x)~y p(x')~zq[x') = 0 

-^(/) + c,p(/) = 0 (3.1.1.) 

-hiy^ + d^qix"^ = 0 

As shown in [2,4] the notions of global asymptotic stability and persistence of the 
system are connected in the sense that the criteria for persistence of the system (i.e. 
the persistence of top predator populations in a simple food chain. Indeed, in our 
system (2.1.1.) each predator acts as a top predator, consequently imply the existence 
of a positive equilibrium. 

Thus, for the existence of positive equilibrium we should have the following conditions: 

(1) Range (g) = c, Range (/?), Range ( h ) = Range (g), and g' 1 and h' both 

exist. 


( 2 ) 


and * 


tf(x) 


p(x) — q(x) are both strictly decreasing functions in the positive cone. 
This for the equilibrium to be globally asymptotically stable (see Remark 3.2.1.). 


(3) -g(0)+c,^(jc)>0 and -h(0)+d ] q(x)>0 where E(x,y,0) and E(x,Q,z) 

are equilibria in the co-ordinate planes (see section 2.2.). These two conditions 
are for persistence of the system (2.1.1) (see [2,3,4]). 
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Remark 3.1.1. 


xf(x) xf(x) xf(x) 

The condition that ^ when 

decreasing also implies global asymptotic stability of the equilibria E(x,y, 0) and 
E(x ,0,z), that is non-existence of limit cycles in the x-y plane andx-z plane. 


(H4) holds 


is strictly 


Thus conditions (1), (2), (3) imply the existence of positive equilibrium. 

We illustrate this by two examples. 

Example 3.1.1. 

Consider the system (2.1.1.) with logistic growth of the prey population and 
Holling type II functional responses. 


x = x(l-x)-y 


2x 
\ + x 


4x 

-z - 

\ + x 


y=y 


~(i +y)+ 


33 2x 
161 + x 


(3.1.2.) 


z = z 


-(l+z)+ 


11 4jc 


9 1 + jt 

Condition (1): is trivially satisfied 
Condition (2): 


d_ 

dx 


p(x) 


= -x < 0, 


d_ 

dx 


*/(*) 

q(x) 


= -i<0, 


0 <x* <x<k 


0 < x* < x < k 
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for condition (3): 

v -Y1 3 

£(*,*,0) = £(j,-,0^ 

£>~,0, ? ) = £(I,0,£ 

- g (0) + ^)=-i + H( p (I)) =3 L>0 

-A(0) + d x q(x) = -\+j q [^j = ^ > °, 

and we have E*(x\y\z*) = £*(.329, .02, .21) a positive equilibrium for (3.12.). 

Example 3.1.2. 

Consider the system: 

x = x(l - x) - yx - 2xz 

y-)[-2(2+y) + 3x] (3.1.3.) 

z = z[-(2+z)+4x] 
condition (2) is satisfied. 

£(*,>>,0) = £^ i,o) 

£(*,0,?M(|a|' 

-h(0) + d x q(x) = 1 > 0, but 
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-#(0) + c iPC*) = -2 + 3.^ = ~ < 0, thus (3) is not satisfied implying that (3.1.3.) 
may have no positive equilibrium. 

No solving 


2 + y = 3x 
2 + z = 4x 



3 


which in not a positive equilibrium of (3.1.3.). 

3.2. Asymptotic Stability of £‘(x*,/,z‘). 

The asymptotic stability of the equilibrium E * (x*, y , z*) for (2.1.1.) is proved 
under the assumption (H4), that is 

q(x) = ap(x). (3.2.1.) 

where a is taken as a positive parameter which will be of significance for the 
Hopf-bifurcation in section 4. 

Theorem 3.2.1. 


[f d (xf(x) 

stable. 

Proof: 


<0, then the equilibrium E*(x\y\z*) is asymptotically 


The entries (a l7 ) of the variational matrix ofthe system (2.1.1.) about 
E*(x\y\z') (as simplification is done using (3.1.1.)) are 
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a " = '(*')* 


xf(x) 


p(x) 

«2I =/c,p(/)>0, 

a 3 , = a</,z*/?(;t*)>0, 


<0, 


a \2 =~p(x')<0, 

°22 =-/$'(/) < 0 , 

a 32 = 0, 


=-ap(x’) <0 

°2i ~ 0 

a 33 = -z*/t'(z*) < 0 


(3.2.2.) 


From (2.2.6.) and (3.2.2.) we have p i >0,p i >0 and p x p 2 -p 3 > 0. 

Thus by the Routh-Hurwitz criterion, Zs*(jc*,y*,z*) is asymptotically stable. 

Remark 3.2.1. 

For the global asymptotic stability of E*(x\y’ ,z'), we define a Lyapunov 
function on {(x,y,z):0< x* <x <k,0<y* <p,0<z* <z<y,P>0,y >0} as 


V{x,y,z)=\(c x +d x ) 




*(*) 

pH 


\dw+ 


It is straightforward to show that 



(3.2.3.) 


xf W 

V(x,y,z) < 0 provided is a strictly decreasing function on q < x* < x < k. 

HOPF-BIFURCATION IN 3-SPECIES SYSTEM 

Hopf-bifurcation is one of the most common ways through which a constant 
steady state loses its stability giving rise to a solution (or a family of solutions) of a 
periodic nature. 

The technique of Hopf-bifurcation relies heavily on the effects induced by the 
parameter incorporated in the system on the stability of the steady state of the system. 
Thus, Hopf-bifurcation is characterized by a stability change of the steady state 
accompanied by the creation of limit cycles. 

The Hopf-bifurcation theorem is a very powerful general tool in establishing 
the existence of a one-parameter family of periodic solutions and in a neighborhood 
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of the bifurcation point, it also gives their periods. However, it does not immediately 
provides any information about the stability of these bifurcating periodic solutions. 

Various successful attempts have been made to ascertain the orbital stability 
of the bifurcating periodic solution such as stability formula as in [6], the procedure 
applied in [5], center manifold theory and the constant "a" value technique in [12], 
and the method of Negrini-Salvadori [7]. Nevertheless, we show by an example that 
the stability of the bifurcating periodic solution seems to be an open question. 

Conditions for Existence of Bifurcating Periodic Solutions 

In this section we obtain conditions in terms of p,, p 2 , p 3 as defined in (2.2.7.) 
for the existence of bifurcating periodic solution emanating from the positive equilibrium 

E*[x\y*,z*) ofthe general 3-species system 

*i =*//(w,>2,•><*), / = 1,2,3 (4.1.1.) 

a being a parameter. 

We shall be seeking conditions on the characteristic equation of the variational 
matrix V(E'), i.e. 

A, 3 + p,A 2 + p 2 A + p 3 = 0 (4.1.2.) 

to have 2 purely imaginary roots and a criterion so as the transversality condition 
holds. 

Suppose (4.1.2.) has one real root A, and two complex conjugate roots: 
X 2 =X r +iX s , A 3 = A r -/A J . 


Thus 


A 3 + p, A 2 + p 2 A + p 3 = (A-A,)(A-A 2 )(A-Aj 

A 3 - A 2 (2A r + A,) +A(|A 2 f +2A r A,)-|A 2 | 2 A,. 


By comparing coefficients of like terms, we have 
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~~P\ — 2X r + X»j 
P 2 + 2X V X , 1 • 
Pi — ”|^ 2 1 ^1 


(4.1.3.) 


Let us consider a, as given in (3.2.1.), to be the Hopf-bifurcation parameter which 
induces (by tunning it about a critical value) the change in stability status of the 

equilibrium E*(x*,y*,z*). 


Now (4.1.2.) has to have two purely imaginary roots if X r = 0. 

Thus from (4.1.3.) we have 

Pi*Pi ~ Pi = ° (4.1.4.) 

Equation (4.1.4.) gives the critical value of a, call it <x c , such that p ] p 2 - p 3 = 0. 
Thus a c satisfies p x ,p 2 -p 3 = 0, 

where p t (i = 1,2,3) denotes the value of p when a is replaced by a f . 

Thus the two purely imaginary roots are ± iw , where w = whit p 2 > 0. 

This is because if X., is a real root of (4.1.2.) then the other two are give by 

^• 2,3 = £ —1 (^i + Pi )^ + Pi) — +^1 Pi ~+~Pi) j *2 


If X , 2 3 are purely imaginary then X,, = -/?, 'and it implies that X . 2 3 = with 

P 2 > 0- 


Next to find a transversality condition in terms of p x ,p 2 ,p 3 such that 




Now eliminating X,, from (4.1.3.) will result in 
Sk\ + 8 X, 2 p, + 2X, r [p , 2 + p 2 ] + p x p 2 -p 3 = 0 


(4.1.5.) 
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Differentiating (4.1.5.) with respect to a, then putting a=a f , and observing that 
X r (a r ) = 0, (P\Pi ~ Pi) a . a = 0, we get 


M + ftL. -ft)..., 


da'' ' ),a ° a < 


A Pi +P2 da 


-T-(P,P,-P >)«_ 


(4.1.6.) 


Hence 


Pi+P2>Q> ^(K( a ))>°(< 0 ) if -^(P iP2-P 3 )<°(< 0 ) (4.1.7.) 

Remark 4.1.1. 

Since for d = — (V( a )) a . a > 0, the eigenvalues cross from left half-plane to 

the right-half plane as a increases, it follows that the equilibrium is asymptotically 
stable for a < a r and unstable for a > a f . 

For d < 0, the eigenvalues cross from the right half-plane to the left half-plane as a 
increases, it follows that the equilibrium point is unstable for a <a c and asymptotically 
stable for a > a f . 

d \ 

This means that for d > 0 i e - ^~(Pi />2 ~ Pi) < if F( a ) ~ PiPi ~ P 3 > then F(a) 

is decreasing at a =a c , which implies that p,p 2 - p 3 < 0 for a > a r , a -a c to be 
small. 

Thus with p, >0 andpj >0, the asymptotic stability of E*(x\y\z*) depends upon 
P\Pi~Pi - Hence, for a <CL e ,p l p 2 -p i > 0 we have £*(/,y*,z*) asymptotically 
stable, while for a>a r , a-a f is small, p,p 2 - p 3 <0, we have E'(x ,y\z*) 
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unstable, whereas for a = a c , p^p 2 - p 2 = 0 there is a Hopf-bifurcation. 

Thus we have proved the following theorem: 

Theorem 4.1.1. 

The interior equilibrium point in the positive cone of the system (2.1.1.) is 
asymptotically stable for a < a c , unstable for a > a c , and there is a Hopf-bifurcation 
at a=tx f ,where a r is given by (4.1.4.). 


However; for —(P\P 2 -Pi) a . a< > 0 [ / - e -'^(^( a )L a( <0 


then theorem 


4.1.1. reads as: 

Theorem 4.1.2. 

The interior equilibrium point in the positive cone of the system (2.1.1.) is 
asymptotically stable for a > a r , unstable for a < a f , and there is a Hopf-bifurcation 
at a = a r , where a f is given by (4.1.4.) 


We illustrate the two theorem by the following examples. 

Example 4.1.1. 

Consider the sysem (The classical Lorenz equations), see [ 6 ] 


x = -ax + ay 


y = -xz + ry-y (4.1.8.) 

z = xy + -bz 

where b and a are positive constants and r to be taken as a Hopf-bifurcation 
parameter. £*(/,/,z‘) = £*(^(r- \),Jb(r - l),(r- 1 )), whitr>l. 


Moreover, p i =]+b+a, p 2 =(r+a)b, p }m = 2ab(r - 1 ). 

P\Pi “£ 3 = 0 implies that the critical value of the parameter r is 
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= ’ J r ( p ' p '~ p '}~r. =^(l + ^ - ^)<0 provided that a >1+6. 


Thus Hopf-bi furcation will occur at r = r c provided that a > 1 +b and this condition 
is the same as in [6]. 

Example 4.1.2. 

Consider the following predator-prey system, see [8]: 

( a k \ 

x = x\ 1- y-z 

V g ; 


y = ^{-y + Pkz) (4.1.9.) 

z = ^(x-z) 


which is a predator-prey model with time lag, where a, k , e, y, P are positive constants 
and a > 0 be taken as the Hopf-bifurcation parameter. The positive equilibrium point 


is E'(x\y,z') = E' 


JL X 

pjt’a* 


1- 


P* 


’ p* 


with p/r - y > 0. 


Thus we get 

Pi = ’ Pi =-^7> A=-T^r(P^-y),hence 
e €p/r e pA: 

P\P 2 ~Pi =^Tp^(«-P^ + Y) = 0 implies a f =pfc-y>0 

( a c is the parameter critical value at which Hopf-bifurcation occurs). 
For the transversality condition we have: 
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-Pi =p^-(p^-y)>0 which implies <0 andhence 

the equilibrium point E'(x\y ,z*) is asymptotically stable for a>a c , unstable for 
n < a c and Hopf-bi furcation occurs at a = a c . 


Stability of Bifurcating Periodic Orbit 

In this section we illustrate by an example that the stability of the bifurcating 
periodic solution seems to be an open question. 

Example 4.2.1. 

Consider a 3-species predator-prey model in which two species (x andy say) 
at the bottom are feeding on each other and at the same time are fed upon by the third 
species ( z ) acting as the top predator of the system. This system can be represented 

as: 


x = jc[l-x + a^-pz] 
y = >'[l-y+ax-Pz] ► 
z = r[-l -z + ax + ay] 


where a, p positive constants (a represents the biomass conversion rate in all cases), 
a is taken as the Hopf-bi furcation parameter. 

The interior positive equilibrium E*(x\y\z *) is given by 


x =y = 


1 + P 


2 = 


3a-l 


1 - a + 2ap 1 - a + 2aP 


with a > 1. 


(4.2.2.) 


Computing p i ,p 2 ,p i , p t p 2 - p 3 , and then from p x p 2 -p 3 = 0, we get. 


a. = -5- with p>2. 
f P-2 M 

Solving the characteristic equation 

A, 3 + p,A 2 + p 2 X + p 3 = 0, and a = a c we get the eigenvalues as 


(4.2.3.) 
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Figure 4.2.1. 


The figure shows the stability status of (4.2.1.) with respect to (is *) as a function 

of the interaction parameters a and p. In the region (a), the positive equilibrium 
point (at which all the three species may coexist) is asymptotically stable. Hopf- 
bifurcating periodic solutions are obtained for all a and p lying on the curve in 
the region (b), the equilibrium is unstable. 

(i) The procedure as in [5] and [9] 

Next we move the origin into the equilibrium is *(x* ,y* ,z*) by co-ordinate 
transformation and performing the parameter transformation (4.2.3.) at the same time, 
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we get the following system which is equivalent to the system (4.2.1.): 




( \ 



P 1 

-1 -P 


f x y 

l 

~2(P-l) 

-(P-2) p ~P(P—2) 

/V 


P-2 

P -1 -p 

P-2 K 

f p 

B-2 B-2 J 


y 

P -{P-2) -P(p-2) 

2P 2p .2 

y 

K Z J 

+ diag(x,y,z) 

y 



P-2 p-2 J 





(4.2.5.) 

The eigenvectors corresponding to the eigenvalues c )j = 1,2,3 are 


S, = co/[l,-l,0], S u = co/[p-2, P-2, 2p]±ico/[tv„(p-lXp-2),w 0 (p-lXP-2),o]. 

Introducing S { together with real and imaginary parts of S i as a new base, we perform 
the co-ordinate transformation 


col[x,y,z] = 


^p-2 

P-2 

2P 


-w, 


.(P-lXP-2) 
% (P-lXP-2) -1 
0 0 


\co, 


l[u,v,w\ 


In the new co-ordinate system (4.2.5.) has the form : 




O 

1 

o 


V 


(p\ > 

V 

= 

o 

o 


V 

+ 

F 2 



o 

o 




F 3 

V J 


where 


(4.2.6.) 


F 1 


F 2 

F 3 


-2w 0 Pmv 

w 0 (P-2)m 2 


/ p 2 (p-2) + 4^ 


P-1 


j«v- w 0 (p-2)v 2 + 


w 


Wb(P-«XP-2) 


~(P - 2)(P + 2)uw + vv 0 (P - 2) 2 vw 


(4.2.7.) 
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The procedure given in [5] and in [9] is as follows; 

The sign of the constant C^a,.) defined in terms of partial derivatives of (4.2.7.) 
determines the stability and direction of the bifurcation of the periodic solution. On 

calculation of C,(a f ) (see Appendix), it is found that Re(C,(a f )) = 0, andhenceno 
information concerning orbital stability of the bifurcating periodic solution. 

(ii) Reduction of the System to the Center Manifold 

We apply (as an alternative method) center manifold theorem which simply says 
that the orbital structure near E'(x\y\z‘) with a replaced by a e is determined by 
the vector field (4.2.1.) restricted to the center manifold* 

The center manifold can locally be represented as 

{(«,v,w) eR i :w = h(u,v) eh( 0,0) = 0,h u (0,0) = A v (0,0)} (4.2.8.) 

we choose h(u,v) as 

h(u, v) = “[A, ,« 2 + 2h n uv + /t 22 v 2 ] + o(3). 

If (w(s), v($), w(s)) is a solution of (4.2.6.) near the origin with a value on the center 
manifold, that is 

w(.y) = h(u(s),v(s )), then 

w(s) -/»„(«, v)ii - h v (u, v)v s 0 (4.2.9.) 

Using (4.2.6.) in (4.2.9.), omitting temis of order at least three and equating coefficients 
to zero, we get 

h u = h ]2 =h 22 -0. This implies w = h(u, v) = 0. 


Thus (4.2.6.) on center manifold becomes 
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u = -w 0 v - 2w 0 Pwv 
v = w 0 u + w 0 (fi-2)u 2 


^P 2 (P-2) + 4 

P-1 


wv-w 0 (p-2)v 2 


(4.2.10.) 


On calculating the constant "a" see [ 12] or the Index "7" (see Appendix), it is found 
that" a " (or I) - 0, and hence no information regarding orbital stability of the bifurcating 
periodic solution. 

(iii) The Negrini-Salvadori Method. 

As our third attempt to study orbital stability of the bifurcating periodic solution, 
we apply the method ofNegijni-Salvadori [7] in which a Lyapunov function is assumed 
in the form: 

V(u,v) = u 1 +v 2 + A(u,v) + B(u,v), 

where A and B are homogeneous polynomials of order three and order four 
respectively. 

A (m, v) = a l u i + a 2 u 2 v + a 2 uv 2 + o 4 v 3 

B(u,v) = b t u 4 +b 2 u i v + b i u 2 v 2 +b 4 uv 4 +b 5 v 4 . 

The coefficient a i and b j are to be determined in such a way that the derivative of V 
with respect to the system (4.2.10.) is 

F(u, v) = G(u 2 + v 2 ) 2 + o(5) (4.2.11.) 

where the constant G is uniquely determined and the orbital stability of the bifurcating 
periodic solution depends on the sign of G. 

It turns out (again) that G- 0, and hence (again) no information regarding 
orbital stability of the bifurcating periodic solution of (4.2.1.) can be deduced. 

APPENDIX 


(I) Calculation of C, (a f ) based on Kazarinoff recipe-summary (see [5,9] Using 
(4.2.7.), i.e. 
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F' = -2w 0 pxy 


F 2 =w 0 ( p-2)x 2 - 


p 2 (P-2) + 4 
P-1 


*y-w 0 (p-2)/ + 


Wo (P-lXP-2) 


F } = -(p - 2)(P + 2 )xz + w 0 (p - 2) 2 yz 

s,,=(f:+f;)/4+i(F>/^)/4=o 


ftl - 



P 2 (P-2) + 4 

2(P -1) 


-2iw 0 


ft. + 2^)/4+!'(^-fJ + 2F' v ) / 4 = pi ^ j- jp 4 2«-.(P ■-1) 

ft, = + Fj, +i(Fl - F;)] [f> + F^]/4+2-'[Fi - F^ +i(f; + £)] 

[Fi-Fj-2tfj]/(4+2<H-.).0 

^l( a r) ~*(i> 2 oi>ll “ ^l^ll | “3 |^ 02 | )/2w 0 +g 2 , 12 - 



^p 2 (P-2) + 4 

2(P -1) 


V 


+ 4w 2 


Hence ReC,(a < .) = 0 

(II) Calculation of the constant "a" 
(see [12, page 277]) 



where all partial derivatives are evaluated at the bifurcation point. 
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Now /' = -2 w$xy, f 2 = w 0 (p - 2)x 2 


P 2 (P-2) + 4 
P-1 


xy-»Jb-2Y. 


Thus a = 0. 
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PERSISTENCE IN A 3-SPECIES 
PREDATOR-PREY MODEL 
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ABSTRACT 

The ecological notion of system persistence in modeling the interaction of two competing 
predatory populations living exclusively on a common prey, is investigated. 

Freedom and Waltman [4], and El-Owaidy and Ammar [2] have discussed the persistence 
of such models based upon the assumption of nonexistence of limit cycles. In this paper, however; 
the nonexistence of limit cycles is proven, first, by global asymptotic stability of equilibria 
through the construction of a suitable Lyapunov function, and second, persistence criteria of 
the system are obtained. 

Key words: Predator-prey system, Lyapunov function, limit cycle, asymptotic global 
stability. 

Mathematics Subject Classification (2000): 92 B05 

INTRODUCTION 

One of the most interesting questions in ecological models concerns the survival 
of all species of the model. This question of persistence becomes more important 
when two or more species compete for a singl e prey speci es. 

The significant concepts of permanence and persistence, both exclude extinction 
of species for all positive initial conditions. In biological terms, persistence means that 
the density of each population remains, asymptotically, above a positive bound 
independent of the initial conditions, i.e.all species stay away from extinction. 

Mathematically, this may be stated in terms of behaviour of solutions of the 
model which represents the biological phenomenon. 

An ecological differential system 
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X, =x i f(x i ,x 2 ,...x ll ) for / = 1,2,...,« 

is said to be permanent or uniformly persistent if there exists a compact set k in the 
interior of R" = {x e R":x i >0 for 1 < i < n such that all orbits end up in k. This 
guarantees that the number of each population x,.(/) is bounded away from zero if 
x ( .(0) > for all/. 

By the term weak persistence, we mean that for all /, lim sup x,(t) > 0 whenever 

oo 

x(0)>0. Under this "lim sup" definition of persistence, a population can frequently 
approach extinction. By the term strong persistence, we mean that for all i, 

Kminf*i(0 > 0 wheneverx(0)>0. 

We have applied the latter definition of persistence as used in Freedman and 
Waltman [4]. This difinition of persistence was reformulated in [4] as follows: 

"A persistence with initial conditions in the positive cone will persist if there are 
no co - limit points of the solution on the boundary of the positive cone". 

This means if □ ( X ) be the orbit through the point X= (x, y, z ) with x>0,y>0, 
z>0, and if Q(X) be the co - limit set of □ (A), then Q(A r ) be interior to the 
positive cone. 

Since, the question of two predator populations competing for a single prey 
population has occupied an important place in ecological literature, so in this paper, 
we have considered a system modeling the interactions of two competing predator 
populations living exclusively on a common prey. For both the predator populations 
we have taken density-dependent death rates and in this sense, wc have extended the 
result of El-Owaidy and Ammar [2]. 

Freedman and Waltman [4] in theorem 4.1., obtained persistence criteria by an 
assumption that for such a system, there are no limit cycles surrounding an interior 
equilibrium in a co-ordinate plane. 

We have actually proved the non-existence of limit cycles by constructing 
Lyapunov functions for the equilibria and have obtained conditions under which the 
equilibria are are globally asymptotically stable. This implies that all the trajectories of 
the system in the positive cone will spiral toward this equilibrium point, that is non- 
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existence of limit cycles. 

Then, we have applied the persistence criteria as given in theorem 4.1. under 
reference, to obtain conditions for persistence of our system. Thus our result of 
persistence seems to be more comprehensive than that in theorem 4.1. of Freedman 
and Waltman [4]. 

The organization of this paper is as follows: 

In section 2, we give the model. Section 3 deals with the global asymptotic 
stability of the problem. Section 4 deals with persistence results and we illustrate 
these by an example. 

THE MODEL AND EQUILIBRIA 

A system modeling the interaction of two competing predator populations y (t) 
and z (/) living exclusively on a common prey x (/). 

x = xf(x)-yp(x)-zq(x) 


y = y[-g{y)+w{xj\ 


( 2 . 1 ) 


z = z[-h(z) + d i q(x)], 

x(0) - x Q > 0,y(0)y 0 > 0,z(0) = z 0 > 0. 


( 

Where 

V 


d_ 

dt 


* c p “i are positive constants, known as food conversion rates. 


We make the following assumptions, which are consistent with models of 
predator-prey systems. 

(HI): /(x); specific growth rate of preyx 

/(0)>0,/ (x) < 0 for all x > 0, there exists k> 0 (the carrying capacity) such 
that f{x) > 0 on 0 < x < k,f ( k ) = 0 and/(x) < 0 on x > k. 

2): p (x); the functional response of the predatory with respect to the prey x and 

p (0) = 0 ,p'{x) > 0 for all x > 0. 
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3): q(x); the functional response of the predator z with respect to the prey x and 

q (0) = 0, q\x) > 0 for all x £ 0. 

(H4): we assume q (x) = a ,p(x), where a is a positive constant. 

5): g (y); density-dependent death rate of the predatory and 

g (y); g'(y)Z0 for all y > 0. 

(H6): h(z)\ density-dependent death rate of the predator z and 
h (0) > 0, A' (z) > 0 for all z > 0. 


Clearly the system (2.1) has equilibiruim point E 0 (0,0,0). By assumption 
(H1), E ( k, 0,0) is also an equilibrium point. We assume that each of the predators 


• / * \ 

y and z can survive on the prey that is there exist equilibrium points E x, y ,0 

V J 


and 


E(x,0,z) such that 



( 2 . 2 ) 


-g\ y + c \P\ * = ° 


and 


x/(x)-zq(x) = 0 

-h(z) + d l q(x) = 0 (2.3) 

w ^ere x,y,x,z > 0 and *x < k,x < k • ( 2 - 4 ) 

GLOBALASYMPTOTICSTABILITYOFEQUILIBRIA and E(x,0,z). 

3.1 Assume (H1)-(H6) hold for the system (2.1) and in a neighbourhood of 


(j.J'.oj i 


in the positive cone, the function 


xf(x) _ 


p( x ) 


is strictly decreasing. 



IDRISH H. ELMABRUK 


149 


Then the equilibrium point is globally asymptotically stable. 


Proof: 

Define a Lyapunov function F(x,y,z ) as [see [1]]: 


x 

V{x,y,z) = j 


1- 


P(*) 

p(w) 


^ f 
+ d, 


*\\ 


1 


«(£) 

?W 


\dw + 


U w-y*^ 


W 


\dw + z 


) 


(3.1) 


Now - 0 an d due to (H2) and (H3), (x, y, z) is positive in the region: 

* * 

0< x < x < k, Ocjxy <p,, 0<z <P 2 , P„ P 2 are positive constants. 


d_ 

clt 


V(x,y,z) = c t 


p{x)-p x 


*/(*) z ji x ) 

p{x) p(x) 


-y 


+rf. 


q{x)-q x 


*/(*) yM 

q{x) q(x) 


-z 


y - yfrgfa)+cA x )] + Z [~K Z )+ d d x )] 


using (2.2), (H4) and with some algebraic manipulations we get 
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+z 


f* 


<xc, p 


*-/>(*) +M°M( z )) 


<0 


Thus £ 



is globally asymptotically stable. 


3.2 Assume (H1 )-(H6) hold for the system (2.1) and in a neighbourhood of (x,0, z ) 

*/(*) 

in the positive cone, the function ^ ^ is strictly decreasing. Then the equilibrium 

point £(*,0,£) is globally asymptotically stable. 

Proof: 

Define Lypunov function (x, y, z ) as: 


V(x,y,z) = j 


1 - 


P(*) 

p(w) 


+d x 


q(w) 


\lw + y + 


2 / * \ 


(3.2) 


Rest of the proof follows as in Lemma 3.1. 
Remark 3.1. 


We consider equilibirium points £ 0 (0,0,0) and £,(£,0,0). 

The eigenvalues of the variational matrix (£ 0 ) of the system (2.1.) about £ 0 (0,0,0) 
are: 


X x = /(0) > 0, X 2 = -g(0) < 0, and X, = -h(0) < 0. (3.3) 

Clearly £ 0 (0,0,0) is a hyperbolic point and is unstable along the x-axis. This implies 
that the prey population x grows near £ 0 . 

The eigenvalues of the variational matrix (£,) of the system (2.1) about 
£,(£,0,0) are: 


X l =kf'(k)< 0, X 2 =-g(0) + c,p(£), X } = -h(0)+d l q(k) 


(3.4) 
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Thus £, (A:,0,0) is asymptotically stable along the x-axis. This implies that the prey 
population * remains in neighbourhood of k. 

Remark 3.2. 


For existence of £^*».y»0j and £(i,0,z) it is necessary that 

g(0) + c x p(k) > 0 and -h{ 0) + d x q[k) > 0 (3.5) 

as it implies increase of predator population y and predator population z. 

PERSISTENCE CRITERIA 

In section 3, we have given necessary conditions for existence of equilibria 

£^x,y,oj an d £(x,0,z), and creiteria for their global asymptotic stability. 

In this section, we shall assume global stability of these equilibria and obtain 
persistence criteria for the system (2.1.). 

first, we prove the following two lemmas. 


Lemma 4.1. The equilibrium j n the interior of the x-y plane is unstable in 

the positive direction orthogonal to x-y plane if 

-/j( 0 )+</JV) >0or - h(0) + d l ap( *1 > 0. 


Proof: 




The proof is immediate upon computing the variational matrix E of system 


V ) 


(2.1) about . 


We have: 
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V 


0 


0 




This if -A(0) + d x q\ x ] > 0 or - h{ 0) + d x ap\ x j > 0. 


we have the result. 

Lemma 4.2. The equilibrium E(x,0,z ) in the interior of the x-z plane is unstable in 
the positive direction orthogonal to x-z plane if -g(0) + c,p(x) > 0. 


Proof : The proof is immediate upon computing the variational matrix v(&) of the 


system (2.1) about £(*,(),£). 

Now, to apply persistence criteria to our system (2.1), we have to check 
hypotheses (Bl)-(B-4) of theorem 4.1. in Freedman and Waltman [4] and 
boundedness of solutions, 
we have: 


G i (x,y,z) = -g(y) + c i p(x) 

G 2 (x,y,z) = -h(z) + d l q(x), 

therefore; condition (B1) is trivially satisfied due to (H1)-(H6). Also notice that 
/?(0) = 0, and p'(x) > 0 

implies p(x) is a strictly increasing positive function. Similarly q(x). 

Condition (B2) is true due to (HI). 

F(0,0,0) = /(0) > 0, F(k, 0,0) -/(*)- 0 
£(xA0)-/-(x)S0. 

Satisfying (B3) there are no equilibria on y-axis or z-axis in y-z plane, for if we 
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suppose that there exists an equilibrium £(0,y, ,z,) in y-z plane which is given by: 
g(y,) = 0 and A(z,) - 0, then this contradicts (H5) and (H6), and satisfying (B4) each 

predator can survive on prey, that is, there exist points and £(jt,0,z) such 

that (2.2) and (2.3) hold. 

Also see remark 3.2. Moreover, we require; 

5): Boundedness of solutions of system (2.1.). We suppose that the functions 

f(x), p(x), q(x), g(y) and h(z) are sufficiently smooth so that the solution 

to the initial value problem (2.1) exists, is unique and continuable for positive 
values of t. Regarding boundedness of solution, see Freedman and Waltman 
[3] and [4]. 

We state and prove the main theorem. 

xf(x) 

Theorem 4.1 Let (B1 )-(B5) hold. Let be a strictly decreasing function and 

-g(0)+c,/?(x) > 0 an d -/r(0)+rf l a^xj>0 (4.1) 

Then the system (2.1.) persists. 

Proof : By (B5) solutions are bounded. 

By remark 3.1 the equilibrium £ 0 (0,0,0) is unstable along the x-axis and 

unstable manifold of £,(£,0,0) is two dimensional. Conditions (4.1) follow from 
lemmas 4.1. and 4.2. Non-existence of limit cycles follows from lemmas 3.1. and 3.2. 

This completes the proof. 

Remark 4.1. 


From(H4), we have 
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a _ V (*) _ ^te of prey consumption per predator z at prey density x 
p(x) rate of prey consumption per predator y at prey density x 

Thus, if we consider a as a parameter then, the system (2.1.) will persit provided 

a> 

Remark 4.2. 

We have discussed persistence criteria for a system modeling the interaction 
of two competing predator populations living exclusively on a common prey. But in 
the same way, the persistence criteria can be obtained for the system modeling 
interactions between two prey populations and one predator population, that is, 

X = xf(x)-zq(x) 

y = yg(y)-zr(y) (4.2.) 

z = z[-h(z)+dtf(x)+d 2 r(yj\ 
and for three-level food web, that is 
x = xf(x)-yp(x)-zq(x) 

y = .vI-gM+ c { p(x)] - zr(y) (4.3) 

z = z[-h(z)+d ] q(x)+d 2 r(y)] 

For construction of Lyapunov functions for the systems (4.2) and (4.3) see [ 1 ]. 
Example 4.1. 

To illustrate the theorem 4.1., consider the system with the Holling type II 
functional response. 
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/, \ 2x 2x 

x = — x)— 3 /--za 


y = y\ 


-(1+^)+ 


1+* 1+X 

33 2x 


16 l+jc 


z = z 


/, x 11 2x 

—(1 + Zj H-Ot- 

v ; 9 \+x 


Here,* = l,^ = ilil 

p(x) 2 


*/M 

Thus vj is a strictly decreasing function for x>0. 
(* ' \ '(1 3 

•H’w 


\ 

= p 

(- 

/ 



p\x 


hi 0) 27 

a > —Vtt- =—. Thus take a = 2 

4; 22 


v 41 2 ' 

£(*,0,z) = £l-,0,- 

To check condition (4.1) 


(4.4) 


-^(o)+^W=-i+^xi=-L>o 
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, n , 11 „ 2 , 44 17 A 

-h{0)+d,ap\ x= -1 +— x2x- = -l +— = — >0 
w 1 U 9 3 27 27 

Theorem 4.1. applies and hence the system (4.4.) is persistent. 
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ABSTRACT 

We extend the well known sequence space l (p) to / (p,v) for o £ 0 and find its K6the- 
Toeplitz dual and continuous dual. Some matrix transformations have also been considered. 


INTRODUCTION 

Let l x and c be the Banach spaces of bounded and convergent sequences 
x = (x k ) with usual norm |x|| = sup|x t |. For p = (p k ) with p k > 0 and sup />*<<» 

k 

Maddox [1] defined the space 



and it was shown that the continuous dual of l(p ), i.e. /'(/>) could be identified 
with l(q ), where (l / p k ) + (1 / q k ) = 1 

Further, in [2] the Kdthe-Toeplitz dual of l{p ) was determined. 

In this paper, we extend the space l(p) to l(p,v>) foro>0 and determine the 
continuous and Kdthe-Toeplitz duals of /(p,u). We also characterize the matrices 
of classes (/(/?, o),/„) and (/(/?, u),c). 
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Definition : 

1.1 For p = (p k ) with p k > 0 and sup p k < 0, we define for u > 0. 

A" -u |x A | Pi <oo 

where l(p) = |.v = (.y*): £ |.v *\ Pk < oo 

and = = ; V*r E ~ q "' lh <00 | 

where £>1 and 1 / p k +1 1 q k = 1 

1.2 Paranormed Space : 

A paranorm g defined on a linear space A", is a function g : X -> R having 
the following properties. 

(i) g(Q) = 0 where Q is the zero element in X. 

(ii) sW = g(--v) for all .r e X. 

(i i i) g(* + y) ^ g( y )+g(j') for al 1 ,v, y e X. 

(iv) The scalar multiplication is continuous, that is A. )( -> X(n —> oo) & 
g(.Y (l -.y) -> 0 as // -» oo for, X n , X e Cand x n , x e X imply that 
g(X„x n - Xx ) -> 0(/i -► oo). 

A paranormed space is a linear space X together with a paranorm g. A total 
paranorm is a paranorm such that g(.r)=0 implies x = 0. 

It is easy to see that l(p,v) is paranormed space by 
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*w-[z*"kr]" 

where M = max(l, sup p k ). 

For; (i) = 

= 0 if and only if x = 0 

(»> *M=[s*'”Kr‘]" 

=[z*>.rf 

-*w 

(iii) M*+^)»[Z*'"k+^r] s 

By minkowski's inequality we have 

^^)s[z*'"krF+[z*-i^rF 

= A(x) + A(>>) 

Hence h(x + y) < h(x) + h(y) 

(iv) IfX„, e CwithA, n ->X, 0 (n->oo) andifx„, a € l(p, u) with 
x„ -> a (n -+ oo) in the sense that h{x n - a) -> 0(n -> oo) 

& h(X„x n - X„a) = £ n~ v \k n .x n - l 0 a\ PH 
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Z"1(^o)k - a )+h(x* -a)+{K -^o)«r 


M 


Now |(X. 0 )| P " < max(l,|X 0 |) = //(say) 

Hence, h(X n x n - X 0 a) -> 0(« -> oo) 

So, /(p,o) is a total paranormed space. 

In the next section we prove that l'(p, u) = M(q,v) i.e. the Kdthe-Toeplitz 
dual of l(p, u) is M(q, u);and /’(/?,o) is isormorphic to M{q,x>). 


CONTINUOUS AND K^THE-TOEPLITZ DUALS 

If (X, h) is a paranormed space with paranorm h then we denoted* the continuous 
dual of X that is the set of all continuous linear functionals on X. If £ is a set of 

complex sequences x = (x t ) then E + will denote the generalized K6the-Toeplitz 
dual of E. 

E* = ja = (a k ): ^a k x k converges for all x e E 

Theorem: 

2.1 Let 1 < p k £ sup p k <oo for all k. Then for u £ 0 /*(/?, u) = M(q, o). 
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Proof: 

Here we use the same analysis as in [2]. Let x e /(/>,o) and a € M(q,u). Therefore 

k k 

< £|a 4 f** u, ‘ / * + ££|jc*r 

k k 

< £| k ' u 

it 

<00 

Hence is convergent. Therefore a e /*(p,u). 

Conversely, let X 0 *** converges and * e /(p,o). This implies that a e M(q,u), 
for otherwise we can determine integers. 

0 = n(0) < n(l) < n(2 ) <. 

such that 

*=»(/-i)+i 

Define a sequence x = (;t A ) by 

x k = (signct k )\a k \ Vi (j + 1)“’* M~/k^ Vk \ forn(j-\)<k<n(j) 

= 0, otherwise. 

Therefore 

iV-hr = 

k-n(j- l)+l 
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, ”U) _,v« 

<(j+\y 2 |y 9 ‘ * (?t 1 Kro+ir p * 

s(; + ir ! 

i.e. x e /(/?,o). But 

«(/) 

S a * jc *= 0' +1 )' 1 
*=»(/-!)+1 

So that Z a *** diverges. Hence a contradictioin and therefore a e A/(</,o). 

* 

Theorem : 

2.2 Let 1 < p k < sup p k < oo . Then /*(/?,u) is isomorphic to M(q,v ). 

Proof: 

It is easy to see that x e (p,u) can be written as x = ^x k e k where e k = (0,0,..., 

k 

0,1,0,...). Therefore for any / e /*(/>,o). 

f( x ) = Z ) = Z **’ where/(^) = a* 

it * 

By theorem 2.1, a e M[q,v) if ^]a k x k converges for x e l(p, o). 

* 

Now, if x e /(/>,o)anda e A/(/>,u); ^a k x k converges and defines a linear 

i 

functional on l(p,v). Furthermore, it is easy to check that 
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whenever g(x) < 1. Therefore ^a k x k defines an element of /*(/>,u). Obviously, 

k 

the map T: l‘(p, u)-> M(q,v) given by T(f) = (a l ,a 2 ) is linear and bijective. 

Hence /*(p,o) is isomorphic to M(q,v). 

MATRIX TRANSFORMATIONS 

Let A = ( a nk ) be an infinite matrix of complex numbers a nk (n,k = 1,2,.) and 

X, Tbe two subsets of the space of complex sequences. We say that the matrix A 
defines a matrix transformation from X into Y, if for every sequence x = (x*) e X 

the sequence Ax = (y„) eY, where Y n = ^a nk x k provided the series is 

convergent. The class of all such matrix transformations from X into Y will be 
denoted by (X, T). 

In this section we characterize the matrices of the classes (/(/?,o),/„) and 
(/(/?, u),c). 

Theorem: 

3.1 Let 1 < p k < sup/?* < oo . Then A e (l(p,v),Q if and only if there exists an 
integer E> 1 such that U n (E)< oo for every n , where 

k 

Proof: 

Necessity : Let A e(l(p,u),Q and* e/(p,o). Therefore <oo for each 

k 

* 

n and for jce/(p,u). Hence (a^)* e l'(p, u) for each n. Since /*(p,u) is 
isomorphic to /*(p,o). Therefore A n (x) = ^ a nk x k is an element of /*(/>,d) . We 
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know that /(/>,o) and l m both are complete and by wilansky [3] that there exists 
a constants K such that 

IK Ml ^ *H 

Hence by the principle of condensation of singularities [4] 

U„(E)< oo, for every n. 

Sufficiency : Let xel(p, o) and the condition U n (E)< oo hold for every n. 
Therefore, for E > 1. 

s X£-*•"“ <«, 

k k k 

Hence A e(l(p,u),l„). 

This completes the proof of the theorem. 

Theorem: 

3.2 Let 1 < p k < sup p k < oo. Then A e(/(p,u),c) if and only if 

(i) for E > 1. U "( E ) < oo for every n \ and 

(ii) a nk —► a k [n -> oo) for every k. 

Proof: 

Necessity : Condition (i) follows as in Theorem 3.1; and condition (ii) follows 
from the fact that e k € l(p, u); 

Sufficiency : Let the conditions (i) and (ii) hold and that x e /(p,u). Therefore 
as in Thbeorem 3.1, we have 
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i.e. ^^„k x k is absolutely convergent for 
* 

each n and for x e l(b, o). Moreover, we can choose an integer k 0 > 1 such that 


( 1 ) 


ft.w= <] 


k ~k {) 


Therefore by theorem 2.2 we have 


(i) Zk* * E T.& Vrt-«*r +1 k«W) 

k'Zk o y 


i/a/ 


where M - max(l,sup /?*). Clearly (1) holds if g* 0 (*) = 0. 
Since 


kzk {i 

Therefore by (1) we have 

!>.<*< -» £ a i*»("“»• “), and hence d e(/(p,u),c). 
This completes the proof of the theorem. 
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ABSTRACT 

In the first part we define the space M'(r) and modified Stielfjes transformation introduced 
by Lavoiene & Misra [1] and Marichev [8] respectively. In the second part we use quasi- 
asymptotic behaviour of distribution to prove initial value Abelian theorem. 

Key Words and Phrases : Modified Stieltjes transformation of distribution, 
quasi-asymptotic behaviour of of distributions. 

AMS Subject Classification Codes (1991): 44 and 46 F. 

INTRODUCTION 

The notion of quasi-asymptotic behaviour of tempered distributions was 
introduced by Zavialov [4], Here we give, first several assertions concerning the 
quasi-asymptotic behaviour at 0 + for elements from S' + . The motivation for our 
investigation made in the first part are theorems of Abelian type for the 
distributional modified Stieltjes transformation, while we shall see in the second 
part that the quasi-asymptotic behaviour at 0 + is more convenient than the ordinary 
asymptotic behaviour of a disbritution especially to prove Abelian type theorems 
for the modified Stieltjes transformation. 


NOTATIONS AND DEFITIONS 

N is the set of natural numbers and N 0 = N<j {0}. R is the Set of real 
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numbers and 0 is the set of complex numbers. SI is the space of 
temperecd distributions with support in the (0,oo). S' is the space of 
distribution of slow growth or tempered distributions which is the dual 
of S. We set R + = {* e R:x > 0} = (0,co) and its closure in R is denoted by 
K ,*:=(-«> ,0]. We denote the partial differential operator by 
D p ,peN 0 . 


Definition 1 [3] : 


A function p : (a,oo)~* R, a eR is called regularly varying function at 
infinity if it is positive, measurable and there exists a real number a such that 
for each x > 0 


lim 

k-+cc 


p(**) 

m 



the number a is called the index of p. 


Definition 2 : 


If function L is regularly varying function with index a = 0 then it is called 
slowly varying function at infinity. Thus L is the positive, measurable on (a,° o) 
from definition 1, we have 


lim 

Jit-* co 


L{kx) 

TW 



Definition 3 (Space M'(r)): 


Let M'(r), r eIR\(~IN) denote the space of all distributions / zSl(IR) 
such that f = t~ r D k F, for some k € INo and locally integrable function F, 
supFc[0,oo) and 

j /F(t) / (/ + (5p /_ * dt < oo for p > 0 

IR 
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Definition 4 : 

The Stieltjes transformation S r (f){s),r eIR\(-IN) is 
function defined by > 

complex valued 


0 ) 

se 0 \(-oo,O], 0 </<oo, reIR\(-IN) 

Modified Stieltjes transformation introduced by Marichev [ 8 ] 

is 

r„(/W)=-=frffi+ £ ) -fWy 

“ vyv " r(a) J A y) y 

0 <y <oo 

It can be written as, 


r(f(x))- 1 } y ^dy . 

■ vw ' it«)l(, + ,r 

( 2 ) 

By putting, r = a - 1 , /(/) = y a '/O') < n (0 
we get, 

By equation (2) & (3) 

(3) 



By interchanging x by z & a by r +1 

r(r+i)r f+1 (/)(2)=5 r (/'/)(z). 

(4) 
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Definiton 5 : Modified Stieltjes transformation 7^ +l 

Let us define the T r+i - transformation r eIR\(-IN). Assume that 

/ e M'(r) and that / = t~ r D k F and J|F(t)|(t + p) _r " l_ * <oo for p > 0, holds. 

« 

Now definition is as follows 

T r+I (f), r eIR\(-N) is complex valued function defined by 

n'+i)u/x*H'+<) . 

o(* + 0 

r elR\(-lN), s €0\(-oo,O], 0</<oo 

Definition 6 : Generalized modified Stieltjes transformation f r+] 

The 7^ +| - transformation of a distribution / eS[(IR) is complex valued function 
T r+i (/) defined by 

r ( r + 1 )^ + i(/)( J )= lim </(')’ Tl(0( 5 + 0" r ' ex P(- vv/ ) > (5) 

W—♦CO 

welR, s e A c r\ e A(s). 

where a is the set of complex number for which this limit exists and /f(.s) is the 
family of all smooth functions, defined on R for which there exists e = e ltJ >0 
such that 0<tj(/)<1, t €/?, rj(f) = l if t belongs to the e - neighbourhood of 
/? + ,r|(f) = 0 if it belongs to the complement of the 2 e - neighbourhood of R + , 
where e >0 is arbitrary if IMS * 0, and 0 < 2e < max Res, if for some 
IMS = 0 and |D p r|(f)| < C„, 

t elR. If r|(f) e /4(.s), s e($\R_) 

(i?_ = (-oo,0]), then r|(f)(s+ 1 )’ r "' exp(-wt) e S(IR) for w eR + , r eR 
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To prove this theorem we use the following definitions, lemma and 
propositions. 

Definition 3.1 ([2]): 

Let A/'(r), r eIR\(-IN ), denote the space of distributions /having 
supports in IR + and admitting decomposition / = t~ r (S + D k F ), where F is a 
function with support in the set F=[a,oo),(cc e/F + ), F(/)f r ~*~‘ is integrable and S 
is distribution having support in IR + YmtB ■ 

Now if / = r r (S + D P F) e AT(r) the modifeid Stieltjes transformation of 
/defined by Lavoine & Misra [1] is 

r(r+i xwx*)=<*('). (»+<r >+1 < 3 *> 


If / e A/'(r 0 ), then we can prove that for all r >r 0 , modified Stieltjes transform 
r(r + l)(f„,/)(*) given by definition of modified Stieltjes equation (5) exists and 


r('- + l)(t 1 /^) = r(r + l)S„ l ( J ) + (r + l)j—) ( 3 .2) 

0 \ s + t) 

For a distribution S e(S[) with support in in a compact the 

modified Stieltjes transform is 




(s+t) r 


(3.3) 


where, \ e C*, supp£ c k yt & £( t ) = 1, t e k e , k* is the e - neighbourhood of k. 
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Definition 3.2 ([2]): 

Distribution / e S[ has the quasi-asymptotic behaviour at o 4 related to 
(l / k) a L( 1 / k ) with the limit g eS' + , g* 0, if for every 4> e S 

!^T(T 7 jfc) < ' k ^ ^ >=< ^ > - 0.4) 


We write in short / l g at 0 + , related to (l / k)° Z(l / *). 

Proposition 3.1 : 

Let /, & / 2 , belong to S' & that/, = f 2 in some neighbourhood of zero. 

^ /, ~gat 0 + , related to ( 1 / k) b L(\/k), then / 2 ~gat 0 + , related to 

(\/k) b L(\/k) too. 

Theorem 3.1 : 

Let f eS' + then, / ~g at0 + related to (1/ kf 1(1/ k) if and only if there 

exists a eR + , a > -b-2, such that F = f a+i */ e L\ oc , F defines a regular 
tempered distribution and 

F (0 ~ Cf b , a+2 {t)L(t), t 0 + . we know that 


Supp Fc/fl t & |F(t)| £(l + ||*||yV , t eR + , a >m. 

Lemma 3.1 : 

Let r - b > 0, b + p +1 > 0. Then 


lim 

*->oo 


*-M) 

1 ( 1 /*) 


oo 



0 


((./nj+ijr* 


Itl*) -+-») 
r( r+ i) 


(3.5) 


L(i ) is slowly varying at 0*, t s IR s s t \ R. ■ 
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Lemma 3.2 : 

Let G e L' l0f (IR), supp G c [0,oo), r > a > - 1 , G(x) - x a as jc -» 0 + and 


I 


c(0 

(*+»r 


W < 00 , z € 


0\(—oo,0] 


Then, 

z-»z'-“r(r+iXr„,c)(z) 

is a bounded analytic function in the domain 6 n B( 0, /?), 0 < e < n / 2, R > 0. 

(Q„ e = {a + Re'*;/? > 0,-n+ €<<f> < 7 t- e},a > 0 & B(0,R) = {z:|z| < R}). 
Theorem 3.2 : Initial value theorem 


Let/ e M'(r 0 ) and that/ ~ g at 0 + related to (l / k) b L( 1 / &) then g = CF i+1 
& for 5 belonging to the closed domain 

Q € = {s €^, arg s e[-7T + t-n - s], (0 < e < n / 4)}. 

We have 


lim (»/*r iVOfe/M) =C EM) 

Z,(l / /r) r(r + l) 


(3.6) 


where, r * 1, r > Z>, and r > r 0 . If L = 1 the convergence is uniform in . 

Proof: 

By definition 3.1, f = t~ r (S + D p F ) where supp/ ci? = [a,oo), (a eIR + ). 
F(x)x' r ~ q ~' is integrable and S is distribution having support in the set 
ke(lR + \B). 
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By proposition 3.1. S-gatO* related to (\ / k) b L(\ / k). 

Hence, S = CD P G , G € L' loc , G ~ CF 6+p+1 as t 0 + (v Theorem 3.1) 

Let k e be the s - neighbourhood of the compact set k and let the function £ be 
from c 00 such that £(0 = 1 . f efc, £(0 = 0 , t ik iE \ |d“£| < C a C(E) 

Using Leibniz Formula, we have, 

S(l) = «')«(') - m>"D r G(,) = D'(|-'5(/)G(r))- (3.7) 

a£p 

a*p 

Where, £ a are functions with their support in k 3a \ /r e ([a,a + 3e]) 

f \ 

v t a -'D k F = D k (t a -'F) - £ D'flJT'F) 

V ; 

Let 5 € 0 , args* 7 t. The modified Stieltjes transformation of the distribution/is 

Itr+lXWX*) = <s(m r%>+(^(by 3.1) 

(s+<) "(i+o 


-(' + 1 ) 


P 


J 


0 


oto 


+ 


a+3e 

I 


a 


g(<)S(<) 




yTUffl 

S i (s +()”'*' 

a *p 




(*+«r 


Now by Lemma 3.1 


lim 

A-+00 


1(1/A:) 


r(r+0(WX*/*)- 


CT(r-b) {r . b) 

I'M) 


(by 3.7) 
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= lim-r—v(r + l) p J 


/r->oo 


L(\/k) 




(38) 


where, 7 /( 5 /£,/),/ e77? + is function with the properties. ^^>0 = 

and 7/(0,/) €l\ 

As r>7> 




(3.9) 


Now we have to prove that the first term of (3.8) tends to zero. 
We have. 


CW-C/„„ 1 (/)i(/) = rf(/)A.,.,(()£(r),< £ /f 


(3.10) 


where d is a locally integrable function and d(t ) -» 0, when t -> 0 + . 

Therefore we can choose & 0 in such a way that for any e > 0,|<7(/ / A:)| < e,/ € k 
and k>k 0 . Then 


**•*> , , n a f C?(r)-/^,(/)!(/) 

*0'*r M ((*/*wr +1 


<7/ 


< 


r^;, . “ f £W Ljfm 

iO '*r '! K(*/*)+»f 


<* 






(3.11) 


We proved that the first term on the right hand side of (3.8) tends to 
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zero when k->ao- 
Equation 3.8 becomes 


! , i:i(iS) r(r+1)(3 ’" ,/)(s/t) ' c 


r M) _ 

r(r + l) 


= 0 


(by (3.9) & (3.11)) 


lim 

k-*°c 


k ~ M) 

L(\/k) 


T(r + \){T r+ J)(s/k) = 


r £M) 

r(r + l) 


-(r-b) 

s 


which is required result, (by result (3.6)) 

For the proof of last statement of theorem. We use lemma 3.2 & Theorem 3.1 
it follows that z r ~ b r(r + \)(T r+i f)(z), is a bounded analytic function in 

Q 0e nB(0,R). Let M(z) = f(r + l)(7^ +l /)(l/z), ze0\(-oo,O]. 


The function w r ~ b M(w), w e f2 0 6 n {w:|w| > 1 / /?}, e > 0, R > 0 is analytic 
and bounded. The first part of theorem 3.2 implies x r ~ b M(x) -> A, as x -> oo, where 


, £H 

r(r+i) • 

This implies that the same holds for w r ~ a Af(w) in the domain 
Q fl6 n{w:|w| > l//?}, a>0, s>0, R>0. 

By the Montel's theorem it follows that w r ~ b M(w ) converges to A uniformly 

in this domain & thus w' ~ b M(w) converges uniformly to A in n o . € as M °°- This 
implies that 

z r ~ b T(r + l)(2^ +1 /)(z) -» A uniformly in as z -» 0. 


Hence the proof of theorem. 
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Lemma 3.3 : 

Let/ e M'(r). Then for x e IR + and s e<£\R_ 


r(r+ixr„,/x«)=^+1)| r(r+2 Xj-„,/X«>* 

X 

Theorem 3.3: 

If f eM'(r) andr(r + 2)(7^ + 2 /)(x)~ .x: -» 0 + , where b<r, 

then r(r + lXr„/X*)~C^*^>I(*), *-+0*. 

\ r ~ b ) 

Proof: 

Let e > 0 be fixed. There exists 8 6 IR + such that for 0 < x < 8 . 

r (r+2Xr„ ! /X*)= c ( 1+ nWK MH iW (3.i2> 

Where e is locally integrable function in the 5 - neighbourhood of zero and such 
that 

|Ce(jc)| <e/ 8 , 0 < x < 8 , with M = C(r - b)~' and N = (r + 1 ). 

We have from Lemma 3.2. 


r(r+iXr„,/Xx) 

NMx‘-'L(x) 



T(r + 2 

Mx b - r L(x) 


4. 

K(n+2Xr„ 1 /X<y/ 


1 Mx‘-L(x) 


Now there exist 8 , 0 < 8 < 8 , such that for 0 < x < 8 


? r(r+2 XWX<>ft 

j Mx‘-'L(x) 


<cii 


Thus, for 0 < x < 8 we have 
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r(r+iXWX*) 

NMx^ r L{x) 



e(tj)Ct~ (r ~ by 'L(t)dt 

Mx b - r L(x) 


-1 + e / 8 


8 t b - r ''L(t)dt 
Mx b ~ r L(x) 


C i 


-l+e/8 


jt b - r -'L(t)dt 

i s^w +E/8 (••i c ^ SE/8 ) 


btx 

C\u b - r -'L(ux)du 


ML(x) 


-l+e/8 


b/x 

ju b - r -'L(ux)du 


ML(x) 


+ e/8 


oo 

ju b - r -'L(ux)du 


C 1 - 


ML(x) 


1+C 


OD 00 

| t b ~ r ~ x L{t)dt | u b ~ r ~' L(ux)du 


Mx b - r L(x) 


+ e/8 


ML(x) 


+ e/8 <e 


. r(r-n)(7;„/X^) rr 

NMx b ~ r L(x) 

We made use of the relation 

oo oo 

J u b ' r ~ ] L(ux)du ~ Z,(*)J u b ~ r ~'du, x -» 0 + 

i i 


(3.13) 


r(r+iX^,/)W 

. NMx b - r L(x) 


(By 3.13) 
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. r(r-MX^,/Xx) 

NMx b ~ r L(x) 

■■■r(r + \){T r+ J)(x)=NMx h ~ r L{x) 

•• r(r+ i)(t; + 1 /)M ~ * -> o + 

(r-b) 

Hence the proof. 


CONCULUSION 

By applying Stieltjes transforms and method given in Pilipovic [3] for quasi- 
asymptotic behaviour of distribution. We obtain "Initial value Abehan theorem for 
Distributional Modifeid Stieltjes Transformations". 
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PEXIDER DIFFERENCE OPERATOR IN L' SPACES 

STEFAN CZERWIK* AND KRZYSZTUF DLUTEK* 

(Received 03.05.03) 

ABSTRACT 

Assume that (G, +) is a complete measurable group with an finite measure X and £ is a 
Banach space. For f g, h : G->E define the Pexider difference operator by P (f, g, h) (x, y) :-f 
(x+y)-g(x)- (h)y, x.yeG. 

In this paper we present some results about such Pexider operator. We prove theorem about 
the Hyers-Ulam-Rassias stability of Pexider operator. Moreover, theorems concerning continuity 
and the existence of continuous inverse, operator to Pexider operator are established. 


INTRODUCTION 

Let £ be a Banach space and let (G,+,Z,X) be a semigroup with zero with a 
complete note identicaly zero measure X and such that X(G) < +oo. 

For any functions f,g,h: G -> E we define Cauchy and Pexider operators 
respectively as follows 

CF(x, y):f(x+y)-f{x)-f(y), 

F{f,g,h)(x,y ): =f(x + y)~ g(x) - h(y). 

Now we prove 

Lemma 1.1 For every f,g,h : G-> E we have 

Cf (x, y ) = P(f ,g,h)(x,y)~ P(f, g,h)( 0, y ) - P(f, g,h)(x,0) - h(0) - g(0), (1) 

Cg(x,y) = P(f,g,hXx,y)-P(f,g > h)(x+y,0) + P(f,g,h)(y,0)-P(f,g,h)(0,y)-g(0) t 

( 2 ) 

* Institute of Mathematics, Silesian University of Technology, Kaszubska 23,44-100 Gliwice, Poland 
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Ch(x,y) = P(f ,g,h\x,y)~ P{f,g,h\0,x+y)+P(f,g,h)(0,x)-P(f,g,h)(x,0)-h(0). 


Proof. We get for x,y € G, 

cf ( '>y) =/(* +y)-f(*)- f(y)=fi x +>')-sW - Ky) +sM+M°) 

-f(x) + h(y) + g(0)-f(y)-g(0)-h(0) 

= P{f >g>h)(x,y) - P(f,g,h)(x, 0) - P(f ,g,h)(0,y) - g( 0) - h( 0). 
Analogously, we obtain 

.v) = g{x + y)- g W ■- g(y) = g{x + y) + h(0)- f(x + y) 

+ f{x + y)- g(*) - h(y) + h(y) + g(0) - f(y) + f(y) - g(y) - h( 0) 

= P{f >g>h\x,y)- P(f,g,h)(x + y,0) + P(f,g,h)(y, 0)- P(f,g,h)(0,y)-g(0). 

The proof of (3) is quite similar. 

Now we establish some denotation and general hypothesis. We define 
v A = {t eG;y + t eA}, A y ={f &G\t+y eA} for ,4 eZ- 
We assume: 

(p) A v eZ and = X.(^) for all A e Z and all A eG, 

(1) y /l eZ andX,[ y y4) = X.(y4) foralM eZ and all y eG. 

The following denotations will be used further on. 

For a complete measure space (A\Z,v) and \<> p< +°o we denote (see [1], 

[2] , [5]) L P V {X,E) := j/:X-*E;/is Z measurable and |f| p < +°oJ, 
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where |f| p := (f x \\f(xf dv(x)) ' for 1 < p < +oo, |/L:=esssup, e J/(f)||. 

Of course, function / and g are equal if they are equal almost everywhere with respect 
to the measure v. 

If f:G->Eandy eG, then we put / v (*) := f(y+x),f y (x) := f(x+y). For 
f:G*G-*E , we denote /°(x) :=/(0,x),/ 0 (jc) :=f(x, 0), f*{x,y) :=f(x+y, 0), 
f + {x,y)-=f{0,x + y). 

We will apply the following 

Lemma 1.2 ([6] Let (G , Z, X) be a complete measure space such that 0 < X(G) < +oo 
and let / e L[(G,E). If ( p ) is satisfied, then for y eG f y eZ|(G,£) and | c / v (/) 

dX(t ) = J G f(t)dk{t ). Moreover, if (/) is satisfied, then for y e Gf y e L[(G,E) and 

STABILITY OF THE PEX1DER EQUATION 

In this section we are going to give the positive answer to the problem of 
Hyers-Ulam-Rassias stability (see also [3], [4]) of the Pexider equation : 

f(x + JV) = g(*) + h(y), x,y eG. 

The problem can be formulated as follows. 

Let f,g,h :G-> E be such that P(f,g,h) eL p ux (G*G,E). Find Constant 

M> 0 and an additive mapping A : G -> E such that the following conditions are 
satisfied 

\f-A\ r SM\\P(f,g,h)\ r , 
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“I Vp 

Note that iHf j\\ p {f>S,h){x,y)\\ P dX(x)dX(y) . 

LcxC 

For the Cauchy equation this problem has been solved in [6] in the case 
X(G) < +oo and in the case X(G) = +oo in [7]. 

By fJ. we understand the completion of the product of the measure X times 

itself. 

The following theorem, which we will apply latter on, can be found in [6]. 

Theorem 2.1 Let G be a semigroup and let (G,L,X) be a complete measure space 
satisfying at least one of the conditions (p ), (/). We assume that the mapping 
s : G x G 9 (*,.y) (x,x + y) e G x G is measurable. Let £ be a Banach space and 

let /:G-+£ be such that Cf eZ£(GxG,£) for a certain l<p<+ oo. If 
0 < X(G) < +oo, then there exists function A :G-+ E such that 

(i) CA = 0 

(ii) f-AeL'{G,E ), 

(iii) \f - A\ p <\(G)-'"\Cf\ p ,. 

(iv) forevery a " :G-+ E satisfying the conditions G4*=0, and f - A* e L{(G,E), 
we have A = A\ 

The next result is the following 

Theorem 2.2 Let the assumptions of Theorem 2.1 be satisfied. Moreover, Let 
f,g,h:G-+E be such that P(f,g,h\P(f,g,h)\P(f,g,h) + eL^G x G,E), 
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P(f,g,h)°,P(f,g,h) 0 eL p } (G,E) forafixed l^/ 7 ^+oo. If 0<X(G)< +°o, then 
there exists A .G-+E such that 

(v) CA= 0 

(vi) f-A,g-A,h-AeL p {G,E ), 

(vii) a) |/-^<A.(G)‘ l/p ||C/|| p , 

b) \ g -A\ p <m i, rgi 

c) \h-A\ i) <X{Gy l ' p \\Ch\\ p , 

(viii) a) if B:G->E, CB = 0,/ - B eL p k (G,E), then B = A, 

b) if D:G->E, CD = 0, g-D &L p k (G,E), then D=A, 

c) if H\G-> E, CH - 0, h- H eL p k (G,E), ihenH = A. 

Proof. At first we quess, in view on Lemma 1.1, that Cf,Cg,Ch e Ifjfi x G,E). 

Therefore, making use of Theorem 2.1, we obtain that there exist A, A,, A 2 :G -» E 
such that 

CA = CA, = CA 2 = 0, (4) 

f-A,g-A„k~A, si?(C,£), 

|*-4|,SX(C)""|ft|„. (5) 


We shall prove that A= A, = A 2 (with respect to the measure X )• 
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Note that we have 

g(x) - A(x) = -P(f, g,h\x,0) - h(0)+f(x) - A(x), 

h(x) - A(x) = -P(f, g,h)(x, 0) - g(0) + f(x) - A(x), 

Whence, using our assumptions, we get 

g-A,h-A eL[{G,E). (6) 

From Theorem 2.1 we have 

if B :G-*E satisfies CB = 0, g-B eL p x (G,E), then# = A,, 

if D .G-+E satisfies CD = 0, h-D eL[(G,E), then Z)=a4 2 , 

Hence, taking B = A and D= A, in view of (4) and (6) we get A - A, = A 2 . 

The conditions (vii) and (viii) we obtain now directly from (5) and (7) 
(substitiuting A = A, = A 2 ). This completes the proof. 

Corollay 2.1 Under The assumption of Theorem 2.2 istead of inequalities (vii) we have 
(ix) a) \f-A\ r i\(G)-' l, \\P(f,g,hl i + \p(f,g,h) , ‘\ r 

b) k- a\ p S +|iv.«.*r|, + 

c> i*-4, s + |p(/.*.*),l, + 
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A p {f^ h W P + Mc) ,/ 1Ko| 

Proof. We will prove the condition a). The inequalities b) and c) one can verify 
analogously. In virtue of Lemma 1.1 we have Cf(x,y) = P(f,g,h)(x,y)- 

P{f>g,h)°{x,y)-P(f,g,h) 0 (x,y) -g( 0) - A(0). Hence we get 

l<y|„ S + M°I, + l*(0l, + 

+ X(G)'''|P(/.g,i)J f + X(G) ! ''[|^0)| + ||A(0f 

Consequely, applying the inequality (vii) a), we obtain the condtion (ix) a), which 
completes the proof. 

PEXIDER LINEAR OPERATOR 

Recall the following theorem contained in [ 6 ]. 

Theorem 3.1 Let G be a non-empty set with a binary operation +. Let be 

a measure space satifying one of the conditions (p) or (/). If 0 < X(G) < +oo, / 

for a certain 1 < p < +oo and the mapping q 2 9 (x,y) -> /(x+y) is measurable in 
the product space, than 

(x) Cf eL'jG xG.E), 

(xi) |/|,SMG)""|C/||,, 

(xii) \Cf\ p £ 3^(G)'"|/| p . 

Now we prove 

Theorem 3.2 Let the assumptions of Theorem 3.1 be satisfied. If f,g, AsZ* fora 
certain l£p£+oo, 0 eG and the mapping Q 1 B (x.^y) f{x + y ), 
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G 2 3 (*» g( x + y), q 2 9 (x, y) -> h{x + y) are measurable in the product space, 
then 

(xiii) Operator: (f,g,h)^P(f,g,h) is linear, 

(xiv) P(f, g,h)L^(G *G,E), 

(xv) If h(0) = g(0) = 0, then 

I/I, + lsl, +W, s + |p(/,g,*)*|[ +|/t/.g.A).| ( ,) + 

3(|f(/.*,A)'| j( +|P(/.g.A),|,) 

(xvi) If A(0) = g(0) = 0, then X(G)‘"'^|f(/'g,*)|, + ||^/.g.*)'|| (i +|A/.g.*),||J 

+|/ > (/.g,A) 0 | < , + h/,g.A)J,S5(|/| i , + |g| i , + |A| p ). 

Proof. If f,g,h,k,l,m e L[, and s eK (set of real or complex number), then we 
get easily 

P(f + k,g + l,h + m)=P(f,g,h)+P(k,l t m), 

P(sf,sg,sh) = sP(f,g,h) 
so P is linear operator. 

Let l r :GxG -> E, r = 1,2,3 be defined by : /,(x,y) := f(x+y), l 2 {x,y) : g 
(x), l 3 (x,y):h(y ). Now we have by the Fubini-Tonelli theorem 

fclj* j\U*,ym*,y)= iH x XM. x *y) s 

G*G G*G 
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and hence l 2 € . By the same arguments we get / 3 € and ||/ 3 ||* = X(G)\h\ p p . 

Similarly in view of lemma 1.2 we obtain 

I/I'=| dl/MT^M = j = M c )"'| aMy)l M x ) - 


X(Cy' - X(C)' , ||/ l ||' i 

GxG 

whence \f\ p = k(G) ,/p ||/ l || p . Therefore P(f,g,/i) = /,-/ 2 -/ 3 el£ and (xiv) has 
been proved. 

Consequently we have also proved that 

+M,+M, - MG)'"(|/| P+ w,+w,> <*> 

Now we shall prove (xv). From Theorem 3.1 we have the inequalities 

I/I, s HG)“"I|C/IU*|, s S M G )'"I CA I,- 

Therefore, by Lemma 1.1, we obtain successively 

w,+w, + w, s +ioi, 


+1 P(f,g,h)-P(f,g,h)‘ - P(f,g.h)° + P(f,g,h) 0 1[ +\p(f,g,h)-p(f,g,h). - 

-P{f.gA + n/.g. A )'lJs 3X(G)-'"[||f(/.*, A | p +|n/.g. A )1 ( , 


+3 


/ , (/.«. A ) 0 |/|/ > (/.«. A )J,} 


Similarly one can easily check that 
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(I/I, +K). h/.S.A)J, S MG)'"(|/I, +W,) 

Adding the above inequalities and the inequality (8), we get (xvi), which proves the 
theorem. 

Let's introduce some definitions. 

DP-.= {(f,g,h) eL'xL'x K-,P(f,g,h)' eL',P(f,g,h\ eL',P{f,g,h)eL' r , 

A(0) = 0,^(0) = 0} 

In DP we define the noim : |/,g.A)|| :=\f\ p +\g\ p +\h\ p ,(f,g,h)eDP, 
and in L' the following one : ||/|| : = >-(G) '''(||F|| f + l /r *||,/) + 

\ F \ + \ F l- F *K- 

Remark 3.1 The set DP is a linear supspace of L P X (G,E) *L p x (G,E)x L P X (G,E). 

Remark 3.2 The following equality DP = {(f,g,h) x L[ x L p x ,(x,y) -> 

f(x + y),(x,y) -> g(x + >»),(/,>^) -> h(x + y) are measurable in the product space} 

is the immadiate consequence of the properties of measurable functions and our 
Theorem 3.2. 

Now we are in the position to present. 

Theorem 3.3 Let the assumptions of Theorem 3.2 be satisfied and g(0) = h( 0) = 0. 

Then the linear operator P :DP P(f,g,h) eL[(GxG,E) is 

continuous and continuously invertible. Moreover, the inverse operator (defined for 
all m e P(DP) has the form : 

P- I m{) = (-X(G)' 1 J c [m(jc,.)-m 0 (x)-m 0 (.)]^(4-X(G)' 1 J 
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-m°(.) + m 0 (.)]^(x) -X(G)' 1 1 c [ro(x,.) - m + (x,.)- m 0 (x) + m°(x)]^(x)) if ( p) 
occurs and 

p-'roQ = (- X ( G r'l g K^)-' m o(-H-A-(G)' 1 J G [m(.,y)-m + 

{•>y)~ m °(y )+^W]^(>) 1 j c [ m ("*y)- m +(-*y) -»*>(•)+^O^M) if 

(/) occurs. 

Proof. The continuity of P follows from Theorem 3.2 point (xvi) and Remarks 3.2. 
The existence of inverse operator P 1 are the consequence of Remark 3.2, Theorem 
3.2 point (xv) and Theorem 3 ,p. 43 from [8]. 

Now we establish the form of p\ Let us assume the condition (p) (in the case 
(1) the proof runs quite similarly). Note that Cf , Cg, Ch e L[(G, E ) and in view of 

Lemma 1.2 we have f(y) = -X(G ) ' J G Cf(x,y)dk(x), g(y) = -k(G) 'j G Cg 

(x,y)dl(x), h(y) = -X(G)“' J G Ch(x,y)dX(x). 

Making use of Lemmal.l we obtain equalities 

f(y )='J G { p (f >S’ h \ x 'y)~ p {f>8’ h ){xfl)- P{f >S,h){ 0 ,y)yk(x), 

g{y) = -X(G)J G {P(f ,g,h\x,y)~ P(f,g,h)(x + y,0) 

+P{f>g,h){yfl) ~ P{f,g,h)(0,y))dX(x), 

h(y) = -X(Gy'j c (P(f,g,h)(x,y)-P(f,g,h)(0,x + y) + 

P(f^g^){0,x) - P(f,g,h)(x,0))dl(x). 

Hence we obtain the form of P 1 . This completes the proof. 
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